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Abstract: This paper presents a worldsheet theory describing holomorphic maps to twistor 
space with N fermionic directions. The theory is anomaly free when J\f = 8. Via the Penrose 
transform, the vertex operators correspond to an J\f = 8 Einstein supergravity multiplet. In 
the first instance, the theory describes gauged supergravity in AdS4. Upon taking the flat 
space, ungauged limit, the complete classical S-matrix is recovered from worldsheet correlation 
functions. 



Contents 



1 Introduction 2 

2 The worldsheet supermanifold 7 

2.1 Automorphisms 9 

2.2 Deformations 10 

3 The twistor string 11 

3.1 Matter fields 11 

3.1.1 The infinity twistor 12 

3.1.2 The action 13 

3.2 BRST transformations 15 

3.2.1 The ghost multiplets 15 

3.2.2 The BRST operator 17 

3.3 Worldsheet anomaly cancellation 19 
3.3.1 Zero modes 21 

3.4 Vertex operators 23 

3.4.1 The N = 8 supergravity multiplet 23 

3.4.2 Picture changing operators 25 

4 Scattering amplitudes in the flat space limit 27 

4.1 A degenerate infinity twistor 27 

4.2 The worldsheet Hodges matrix 29 
4.2.1 Self-dual M = 8 supergravity 33 

4.3 The conjugate Hodges matrix 34 

4.4 The tree-level S-matrix 37 

5 Discussion 38 

5.1 The SL(2; C) system 38 

5.2 Higher genus 39 

5.3 Boundary correlation functions in AdS4 40 

5.4 Other issues 42 

A Some properties of algebraic /37-systems 44 



- 1 - 



1 Introduction 



The most influential scattering amplitude in Yang-Mills theory is undoubtedly the Parke- 
Taylor amplitude [1] 

^°-(l,2)(2,3)-..(u,l) • (L1) 
It describes the tree-level (colour-ordered) scattering of two gluons i and j of negative helicity 
and n — 2 gluons of positive helicity, each of momentum = AjAj. Its talismanic status rests 
on two pillars. First, its extraordinary simplicity assures that scattering amplitudes are far 
more managable objects than could be expected from momentum space Feynman diagrams, 
encouraging us that their structures and properties will repay our close attention. Second, 
it provides a remarkably fertile base for deeper explorations of the full Yang-Mills S-matrix. 
In it, one already sees hints of the twistor action for Yang-Mills [2, 3] and the associated 
MHV diagram formalism [4-7], the germ of the Grassmannian formulation of all on-shell 
diagrams [8-10], and the amplitude's factorization properties — a crucial ingredient of BCFW 
recursion [11] — laid bare. Finally, and of particular relevance to the present paper, (1.1) is 
an avatar of Witten's representation [12] 1 

r d 4(k+2)|4(k+2) z 1 n 

= J vol(GL(2;C)) (12)(23)-.(nl) U^™) ^ 

of the n-particle N k MHV amplitude in TV = 4 SYM as an integral over the space of degree 
k + 1 rational curves in twistor space. Witten obtained this form by generalizing Nair's 
interpretation [13] of the Parke- Taylor amplitude in terms of (the leading trace part of) a 
current correlator supported on a holomorphic twistor line. 

This paper is concerned not with Yang-Mills theory, but with gravity. An expression for 
all n-particle tree-level MHV amplitudes in gravity was found by Berends, Giele &; Kuijf [14] 
only two years after the discovery of the Parke- Taylor amplitude. Despite this, the gravita- 
tional S-matrix has proved more resistant to study than its Yang-Mills counterpart. Although 
gravity amplitudes admit a BCFW expansion [15-18], actually carrying it out leads to ex- 
pressions for n-particle N k MHV amplitudes whose structure is as yet unclear [19]. Applying 
Risager's procedure to the BGK amplitude [20, 21] leads to an MHV diagram formulation 
that fails when n > 12 [22, 23], while the twistor action for gravity tentatively proposed 
in [24] does not appear to extend naturally to J\f = 8 supergravity. Considering that M = 8 
supergravity is supposed to be the simplest quantum field theory, it has been remarkably 
difficult to grapple with. 



throughout this paper, (ij) denotes the SL(2; (C)-invariant inner product eapof-ay of the homogeneous 
coordinates cr— on an abstract curve E of genus zero. Z denotes a holomorphic map Z : E — > PT to A/"-extended 
supertwistor space, here with N = 4. In (1.2) this map has degree k+ 1. Ai(Z) are twistor representatives of 
the external wavefunctions. See [121 for further details. 
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The situation changed dramatically with Hodges' two papers [25, 26]. Hodges showed 
that the tree-level MHV amplitude for gravity could be reformulated as 

M n ,o = (iJ) 8 det'(R)5 4 (^2 P }j , (1.3) 

where H is the n x n symmetric matrix with entries 

H« = ]H ^ = -£MHM. (1.4) 

where |a) and |i>) are arbitrary spinors. The diagonal entries are nothing but the charac- 
teristic gravitational 'soft factors' for the i th particle [14, 27, 28]. H has rank n — 3, and 
det (H) is any (n — 3) x (n — 3) minor of H, divided by the permutation symmetric combi- 
nation (j"i,r2)(r2,r3)(r3,ri) corresponding to the removed rows and also by a similar factor 
( c i) c 2)(c2, C3)(c3, c\) for the removed columns. Hodges' representation has many remark- 
able properties. Chief among these is that Bose symmetry in the external states is achieved 
through determinant identities rather than through an explicit sum over permutations. 

Like the Parke- Taylor amplitude (1.1), (1-3) provides an inspirational starting point from 
which to launch deeper investigations of the gravitational S-matrix. It opens up a path by 
which to approach gravity on its own terms. In particular, unlike the BGK form, (1.3) makes 
no mention of any cyclic ordering that is an artifact of trying to fit gravitational pegs into 
a Yang-Mills hole. See [28, 29] for deconstructed forms of the Hodges amplitude that were 
known previously, and [30, 31] for a graph-theoretic explanation of the relation between them. 

One outcome of these investigations was given in [32], where it was conjectured that 
arbitrary n-particle N k MHV tree- level amplitudes in M = 8 supergravity could be represented 
as 

r d 4(k+2)|8(k+2) Z n 
Mn ' k = J vol(GL(2;C)) det ( H ) det ( H ) n^(^)H^) . (1.5) 

This form was obtained by interpreting (1.3) in terms degree 1 holomorphic maps from a 
Riemann sphere E into twistor space, and then generalizing to higher degree maps. Thus it 
bears the same relation to (1.3) for j\f = 8 supergravity as (1.2) does to (1.1) for j\f = 4 SYM. 
In (1.5), H is the n x n matrix of operators 



1 

W) 



d_ _d_ 

dm ' dfij 



k+2 



that act on the external wavefunctions hi(Z), generalizing (1.4). Here, the a r are any k + 2 
points on E. The factor of det'(H) in (1.5) is any (n — k — 3) x (n — k — 3) minor of H, divided 
by the Vandermonde determinant 

l°n ••• <7 nc+sl = II ( r * r i) ( L7 ) 

i<j 

*jS{ removec i} 
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of the worldsheet coordinates corresponding to the removed rows, and a similar factor for the 
removed columns. The combination det (H) is independent of the choice of minor. We shall 
call HI 'the worldsheet Hodges matrix', or often just 'the Hodges matrix'. Similarly, H v is the 
n x n matrix with elements 

TnrV _ (A(o-/),A((J m )) TnrV _ (H^l) ' H a rn)} n TT 2 ( a s m ) TT ( k 



Im 



fori#m , K= - ^-~ n n 



(1.8) 

where the a s are any n — k — 2 points. The factor of det'(H v ) in (1.5) is any (k+ 1) x (k+ 1) 
minor of H v , divided by the Vandermonde determinant 

K •■•ffi 1e+x | = ]J ( r l r m) (1-9) 

i,m£{remain} 

corresponding to the rows that remain in this minor, and again by a similar factor for the 
remaining columns. Again, though it is not obvious from our current description, det'(H v ) 
is completely permutation symmetric in all u sites. Under a parity transformation of the 
amplitude, det'(Hf) and det'(IHI v ) are exchanged [33-35] (up to a Vandermonde factor involving 
all n points), hence we shall call H v the 'conjugate Hodges matrix'. When k = 0, the 
Vandermonde determinant (1.9) should be taken to be unity and det'(H v ) itself is almost 
trivial. This is why it is invisible in (1.3). 

The conjecture that (1.5) correctly describes the complete classical S-matrix of N = 8 
supergravity was proved (to a physicist's level of rigour) in [33], where it was shown that (1.5) 
obeys the twistor space form of BCFW recursion [36-39] at the heart of which is the require- 
ment that the amplitude has the correct behaviour in all factorization channels. (1.5) has 
also been shown to possess the correct soft limits [35]. For preliminary investigations of a 
Grassmannian representation of (1.5), see [33, 34]. Using this Grassmannian at k = 0, an 
investigation of the MHV diagram formalism for gravity very recently been carried out in [40] ; 
excitingly, it has the potential to overcome the limitations of Risager's method. A different 
(presumably equivalent) generalization of Hodges' form to higher degree maps can be found 
in [41, 42]. 

The most striking property of the representation (1.5) is that it exists at all. The unfath- 
omable morass of Feynman diagrams that contribute to an n-particle gravitational scattering 
process miraculously conspires to ensure that the tree amplitude lives on a rational curve in 
twistor space! At MHV, this fact was originally seen by Witten in [12] using the BGK form 
of the amplitude 2 . It was also shown to hold for the 5-particle NMHV amplitude in [43] 3 . 
The existence of (1.5) means, first and foremost, that all gravitational tree amplitudes live 
on higher degree rational curves in twistor space. 



2 The derivatives in the Hodges matrix (4.2) are responsible for what was called 'derivative of a 5-function 
support' in [12]. 

3 The non-trivial statement here is just that the 5-point NMHV amplitude has support on some CP 2 C PT. 
Any five points on a plane define a conic. 
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Why should the gravitational S-matrix know about these curves? The answer pursued 
here — really, the only conceivable answer — is that four-dimensional gravity is a twistor 
string theory. The purpose of the current paper is to find this twistor string. Specifically, 
over the course of this paper we shall construct a worldsheet theory that localizes on holo- 
morphic maps to N = 8 supertwistor space and whose vertex operators correspond via the 
Penrose transform to a linearized N = 8 supergravity multiplet. We shall obtain (1.5) from 
worldsheet correlation functions of this theory at g = 0. This work thus provides the theo- 
retical framework in which (1.5) should be understood. See [31, 44-51] for earlier attempts 
to understand Einstein gravity in the context of twistor strings. 

It is clear that (1.5) possesses a rich and intricate structure. What clues are there to guide 
us to the underlying theory? The main prompt follows from a simple counting that also helped 
stimulate the discovery of (1.5). At g loops, n particle gravitational scattering amplitudes are 
proportional to the K n+2g_2 , where k is the square root of the Newton constant Gn and has 
dimensions of (mass) -1 . Since the amplitude itself is dimensionless, these dimensions must be 
balanced by kinematic factors. But on twistor space, the only object that fixes a mass scale is 
the infinity twistor — an antisymmetric bitwistor whose presence breaks conformal invariance. 
With flat space-time, the infinity twistor appears in two separate guises, corresponding to the 
( , ) and [ , ] brackets familiar from spinor momenta. With n particles at g loops, the twistor 
space amplitude needs to contain precisely n + 2g — 2 factors of ( , ) and [ , ] in total. Under a 
parity transformation ( , ) and [ , ] are exchanged, along with the numbers n± of positive and 
negative helicity gravitons participating in the scattering process (in the pure gravity sector). 
We deduce that the twistor space amplitude must be proportional to 

n + + g — 1 = n — k+q — 3 factors of [ , 1 and 

(1.10) 

n_+g— l = k+ g + l factors of ( , ) . 

Note that the symmetric choice (n + 2g — 2)/2 is not possible since u may be odd and 
the integrand is rational. The fact that n+ and n_ respectively go with [ , ] and ( , ) is a 
convention fixed by comparison with (1.3). This dependence is certainly realized in (1.5), 
where det (H) is easily seen to be a monomial of degree n — k — 3 in [ , ] whereas det'(IHI v ) is a 
monomial of degree k + 1 in ( , ). At higher loops, (1.10) is compatible with all factorization 
channels of g-loop N k MHV amplitudes. 

The key question is to ask what sort of worldsheet objects could be responsible for this 
behaviour. For [ , ], the answer we find is that k — g + 3 of the vertex operators are fixed and 
do not involve the infinity twistor, whereas the remaining n — k + g — 3 are integrated and 
are linear in [ , ]. To achieve this, the worldsheet theory will involve a field with (generically) 
k — g + 3 zero modes whose fixing is associated to the Vandermonde factor (1.7) in det'(lHI). 
The dependence on [ , ] in the integrated operators is introduced by the BRST operator 
responsible for the descent procedure. The infinity twistor in the form [ , ] really endows 
twistor space with a twisted holomorphic Poisson structure, and the BRST operator we 
arrive at is somewhat reminiscent of those found in Poisson sigma models [52, 53]. 



- 5 - 



The theory also contains a different field with (generically) k + g + f zero modes whose 
fixing provides the Vandermonde determinant (1.9) in det'(H v ). The vertex operators associ- 
ated to this field have a rather different character that may be motivated as follows. Firstly, 
notice that the previous integrated and unintegrated vertex operators have apparently al- 
ready used up all the available n insertion points. (1.5) allows us to choose the minors of 
EI and H v completely independently, so there seems to be no compelling reason to place the 
'additional' operators preferentially with either type of insertion mentioned in the previous 
paragraph. Secondly, although det'(H v ) knows about the holomorphic map Z : £ — > PT, it 
is completely oblivious to the external world. Unlike det'(H) which knows about the external 
states through the derivative operators in (1.6), nothing in the definition of det'(IHI v ) cares 
what we choose for the wavefunctions hi(Z), nor even how many particles are being scattered. 
All this is strongly reminiscent of 'picture changing operators' of the RNS superstring. See 
e.g. [54, 55] for an introduction to these operators. We shall indeed find picture changing 
operators in our theory, and inserting k + g + 1 of these leads to the requisite dependence on 

(,)• 

The characterization sketched above may sound worryingly piecemeal. On the contrary, 
the whole theory flows naturally from a single structure: the worldsheet carries a certain 
exotic twisted supersymmetry. All the required objects fit into geometrically meaningful 
worldsheet supermultiplets — properly understood, the theory contains only three different 
fields. The action and BRST operator £11*6 clS simple as one could wish. 

The outline of the paper is as follows. In section 2 we describe the worldsheet superman- 
ifold whose fermionic symmetries and moduli lie at the heart of the whole construction. The 
actual worldsheet theory is a relative of Berkovits' formulation [56] of the original twistor 
string, and is described in section 3. (Some readers may prefer to begin with this section.) 
Here we begin with a description of the worldsheet fields in section 3.1 and BRST operator 
in section 3.2. After a brief diversion, we proceed to show in section 3.3 that the model is 
anomaly free if and only if the target twistor space has N = 8 supersymmetry. We conclude 
our discussion of the general theory in section 3.4, presenting the vertex operators of the 
model and explaining their relation to an N = 8 supergravity multiplet. Section 4 contains 
the derivation of the complete flat space S-matrix of classical M = 8 supergravity (1.5) from 
correlation functions of vertex operators on the worldsheet. The model of section 3 describes 
SO(8) gauged supergravity on an AdS4 background in the first instance. Thus, before em- 
barking on the S-matrix calculation, in section 4.1 we show in section 4.1 how to rescale the 
fields so as to take the flat space limit. The Hodges matrix (1.6) and the conjugate Hodges 
matrix (1.8) have a different origin on the worldsheet. They are obtained in sections 4.2 
and 4.3, respectively. It is worth pointing out immediately that the present model is more 
successful than the original twistor strings [12, 56] were (as a theory of pure J\f = 4 SYM) in 
at least one respect: the worldsheet correlator we consider leads inexorably to (1.5) and only 
to (1.5). No terms are ignored or discarded by hand. Our work suggests many promising 
avenues for future research. We conclude in section 5 with a brief discussion of some of these. 
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2 The worldsheet supermanifold 



In this section we describe the geometry of the worldsheet supermanifold X on which the 
twistor string theory is based. See e.g. [57-59] for much more information on complex super- 
manifolds. 

Let £ be a closed, compact Riemann surface of genus g. We extend £ to a complex 
supermanifold X of dimension 1|2 by picking 4 a line bundle C — > £ of degree d > and a 

1/2 

choice of spin bundle . We then define X to be the split supermanifold whose tangent 
bundle TX is 

TX = T£©T>, (2.1) 
where T> is the rank 2 fermionic bundle 5 



v^n(c 2 ®K y 1/2 ® C) (2.2) 



over £. We will often say that objects taking values in Kj-. <g> C q have spin p and charge q. 
Thus, sections of T> have spin — | and charge +1. 

For a local description of X, we cover the bosonic Riemann surface £ by open coordinate 
patches U a and let {U a } be the corresponding cover of X. Each such U a is (an open subset 
of) C 1 ' 2 , so we may describe X locally in terms of one bosonic and two fermionic holomorphic 
coordinates z\9 a , with a = 1,2. The fact that X is a split supermanifold means that on 
overlaps the coordinate transformations are 



{g a p{z p )) a b 9% 



v2-3) 



where the transition functions f a p and g a p on U a n Up depend only on the bosonic coordinate 
z, not (z\8). To identify these transition functions, suppose we write an arbitrary section 
V : X — > TX of the tangent bundle of X as 

V = V z (z\0)§- z +V a (z\9)-^. (2.4) 

Recalling that TX = TS © T>, we see that V z d z is a section of TS (written in terms of the 
local basis d/dz) whilst V a d a is a section of T> (written in terms of the local basis d/d6 a ). In 
order to compensate the transformations of the basis 

i = 7^ and 4 =( ^ v 4 (2,5) 

that follow from (2.3), the components V z and V a must transform as 

V z (z a \e a ) = f , a pV s (z fi \6p) and V a (z a \9 a ) = {g a p)\ V\zp\6p) (2.6) 

4 Eventually, the twistor string path integral will include a sum (or integral) over all such choices. 
5 Here, 77 is the 'parity reversing functor' whose role is simply to remind us that T> is fermionic. 
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on overlaps. But since V C°l 2 © T^ /2 © C, we have 

GmA = X (fc^ft (2.7) 

where each component of the 2x2 matrix h is a transition function for sections of C. 

Because X is a split supermanifold, it can be viewed as the total space of a bundle over 
£ — indeed, this is just what is meant by the transformation laws (2.3). To identify this 
bundle, note that by (2.7) and (2.6), the coordinates 9 themselves transform as components 
of a section of T>. But since the coordinates on a bundle transform oppositely to the bundle 
itself, we find that X is the total space of T) y — > S, where P v is the dual of T>. To say this 
differently, functions on X are superfields <&(z\9) that may be expanded in the usual way as 

*(z|0) = <f>(z) + 9 a ^ a {z) + l -e ab 9 a 9 h i{z) (2.8) 

where <p(z) is a function on S, the ip a are a pair of functions on S with values in © C~ 1 

and of opposite Grassmann parity to (f>, and £ is a function on £ with values in Ky, © C~ 2 . 
Thus, the structure sheaf of X is Ox = Oe(A*X> v ). More generally, if $^(^10) is a section 
of i£g (8) C q , then it may be expanded in terms of fields cj)^ >,q \ ip^'^ and ^ p ' q ^ on S, of (spin, 
charge) = (p, q), (p + |, g — 1) and (p + 1, g — 2), respectively. 

To give an example that will be important later, suppose that £ is the Riemann sphere. 
On CP 1 , the bundles K~ l l 2 and C are uniquely determined to be O(l) and 0(d), respectively. 
Thus V = C l 2 ©O(d+l) and the supermanifold X may be identified as the weighted projective 
superspace WCP( 1]1 i ( j +1j( j +1 ) with homogeneous coordinates (<7-|$ a ) obeying the scaling 

(a^ a ) = (ra^\r d+1 $ a ) for any r £ C* . (2.9) 

In this function G Ox may be expanded as 

$(<x|tf) = </>(a) + + ^e ab r^(^) , (2.10) 

where 0, ^ a and £ have homogeneities 0, — (d + 1) and — 2(d + 1), respectively under the 
scaling (2.9). 

Returning to the general case, the cotangent bundle T V X to the supermanifold is just 
the direct sum T V X = K s © V y dual to (2.1). Thus the holomorphic Berezinian Ber(X) of 
X is 

Ber(X) = Ber(Xs) © Ber(£> v ) . (2.11) 

To compute this, recall that for an even parity (bosonic) bundle Ber(S) = Det(-B) — the top 
exterior power of B. However, for an odd parity (fermionic) bundle Ber(i ? ) = Det(ilF) -1 , 
where the power of —1 appears as a consequence of the fact that in Berezin integration 
the integral form d#id#2 transforms oppositely to the differential form d£i A d^2 involving 
variables Q a that have same quantum numbers, but opposite Grassmann parity to 9 a . (In one 
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dimension, this is just the familiar statement that since J d9 = 1 by Berezin integration, if 
9 — > g(z)8, we require that the integral form d9 — > g~ 1 (z)d6.) Thus we have 

Ber(X) = Ky, <8> Det(77'D v ) _1 = <8) <8) £~ 2 ) 1 
^ £ 2 , 

where we used the definition (2.2) in the second step. 

When we come to write the worldsheet action in section 3, we will need a top holomorphic 
integral form on X. Since the Berezinian of X is isomorphic to £ 2 , the product Ber(X) ® C~ 2 
is trivial. Thus it admits a global holomorphic section that we write as d 1 ' 2 ^. For example, at 
genus zero d 1 ' 2 ^ = (o"d<7)d$idi92 in terms of the homogeneous coordinates (a— |$ a ) introduced 
above. We can treat d 1 ' 2 ^ as a top holomorphic (integral) form on X of charge —2. In order 
to construct an action, this charge must be balanced by the worldsheet Lagrangian L, so that 
d 1)[2 zL(z\6) may be integrated over X. 

Let us close this subsection with a couple of remarks. As usual for complex supermanifolds 
(and as on the twistor target space CIP 3 '"^) we take X to be a cs manifold [57, 58], in the sense 
that the antiholomorphic tangent bundle is TX = T£ and so has rank 1|0. Antiholomorphic 
fermionic directions simply do not exist — all operations with the fermions will be purely 
algebraic. Finally, we note that X is not an M = 2 super Riemann surface (see e.g. [59]), 
because our choice of TX means the distribution T> is integrable in the sense that {T>, D} C D. 
Indeed, the usual superderivatives 

on an M = 2 super Riemann surface do not make any sense for us, because the second term 
in each expression has different charge from the first and hence is forbidden. Exactly these 
forbidden terms are responsible for the non-integrability of the odd distribution on an N = 2 
super Riemann surface. 

2.1 Automorphisms 

We now consider the symmetries of X as a complex supermanifold. On a local patch U a , 
as usual these are generated by vector fields V £ 0,°(U a , TX|^ ). We will be particularly 
interested in the symmetries of the distribution T> — these are generated by the vector fields 
V € f2°([/ Q , T>\q ) that act trivially on S. (This restriction would not make sense on a super 
Riemann surface, but does make sense on X precisely because V is integrable.) 
From (2.4) we can write 

Va ^W- = (v a (z)+O b R a b (z) + l 2 e bc e b 6^{z)^ A (2.14) 

where the components v, R and v have (spin, charge) = (— |, 1), (0, 0) and (+5, —1), respec- 
tively. Because the fermionic distribution T> is integrable, the anticommutator {Vi, V2} of any 
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two such vector fields again lies in D. A short calculation shows that the component fields 
obey the algebra 



where 



[vi , R2] = V12 [vi , R 2 ] = vu 

[Ri,Rz] = R\2 {V1,V2} = R12 , 

v? 2 = {R2TA v a 12 = (R2) a b v\ - tr(7? 2 )v? 

(RwTb = (R2R1 - RlR2Tb (R'nTb = ~ v 2 V U ■ 



(2.15) 



(2.16) 



with vj, = efe c v c . All other commutators are zero — in particular, {vi, V2} = since it is of 
order (#) 3 which must vanish. 

If we decompose the gl(2; C) matrix R as 

R\ = \8\r + r\ (2.17) 

where the traceless, symmetric matrix r a b takes values in sl(2; C) while r = tr(i?) takes values 
in gl(l ; C). R may be interpreted as generating a gauge transform of C 2 (g>£ so that r generates 
gauge transformations associated to the determinant C 2 . Equations (2.15)-(2.16) then reflect 
the fact that the v a transform in the fundamental representation of SL(2; C) and have charge 
+1 under C, whereas the v a transform in the fundamental of SL(2;C) but have charge — 1 
under C. 

Later, in section 3.2 we shall introduce a (bosonic) ghost multiplet in 77il°(X,D) cor- 
responding to (2.14). This algebra will then be interpreted as the gauge algebra of our 
worldsheet theory. Zero modes of the ghost multiplet live in H°(X,T>), parity reversed, and 
correspond to globally defined ^-preserving infinitesimal automorphisms of X as a complex 
supermanifold. 

2.2 Deformations 

As for a usual complex manifold, infinitesimal deformations of X as a complex supermanifold 
are parametrized by elements of H 1 (X.,T'K). This cohomology group is the holomorphic 
tangent space to the moduli space of X. Again, we will be interested in the moduli of X 
associated to deforming the choice of distribution T> whilst leaving S fixed. Infinitesimally, 
these are described by iJ 1 (X,P) C i? 1 (X, TX). On a supermanifold, the dualizing sheaf is 
the holomorphic Berezinian, so the deformations are Serre dual 6 to H°(X, Ber(X) <g> T> v ). We 
showed in (2.12) that Ber(X) = C 2 , so using (2.2) we can identify this cohomology group as 

H°(X, Ber(X) <g> V y ) ^ 77F°(X, C 2 ® K^ 2 ® C) , (2.18) 

where again the symbol 77 denotes that the fibres are fermionic. This group is non-trivial 
provided only deg(£) > 0, so X has odd moduli even at g = 0. In section 3.2 we shall introduce 
a (bosonic) antighost multiplet valued in 17° (X, C 2 ® K^/ 2 (g> C). Zero modes of this antighost 
live in (2.18), parity reversed, and so by Serre duality can be paired with deformations of 



See e.g. [59-61] for a discussion of Serre duality for complex supermanifolds. 
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the odd moduli. This is the usual mechanism by which the RNS superstring provides a top 
holomorphic integral form on odd moduli space; see e.g. [55, 59]. 

3 The twistor string 

In this section we define the worldsheet theory that will provide a twistor description of 
Einstein supergravity. After introducing the fields and explaining their geometric meaning, 
we study the gauge and BRST transformations naturally associated to the structure of the 
worldsheet supermanifold X. The model is chiral, and we show that all (local) worldsheet 
anomalies vanish if and only if the target space has W = 8 supersymmetry. We then con- 
struct vertex operators in the BRST cohomology, finding that they correspond to an M = 8 
supergravity multiplet. We have just seen that X has odd moduli even at genus zero. We 
construct the associated 'picture changing' operators. 

3.1 Matter fields 

To define the worldsheet model, we first introduce four bosonic and AT fermionic fields Z 1 
(where I = 1, . . . , 4|1, . . . ,Af). Each of these are scalars on X of charge +1. In other words, 

2eS]°(X,C 4|Af ®£) (3.1) 

where C is the same degree d line bundle used in the definition (2.2) of P. In the first instance, 
Z defines a smooth map Z : X -> C 4 ! 4 . The twisting by £ means that this map is defined only 
up to an overall non-zero complex rescaling, so Z really defines a map Z : E — > CP 3 '"^. The 
Z 1 then represent the pullbacks to X of homogeneous coordinates on this projective space. 

Saying that Z is a map from X, rather than from E, simply means that it is a worldsheet 
superfield. As in (2.8), we define its component expansion to be 

z\z,e) = z\z) + e a P i(z) + \<ro a Y\z) (3.2) 

in terms of fields (Z 1 , p^Y 1 ) on S. Since each 6 a is a fermionic coordinate of (spin, charge) 
= (—2,1) we see that Z 1 , like Z 1 , is a scalar on E of charge +1, p l a are a pair of uncharged 
(|, 0)-forms, while Y 1 is a (l,0)-form of charge —1. If the index I denotes a bosonic direction 
in CP 3 ' - ^, then Z 1 and Y l are bosons while p^ are fermions. This is reversed when / denotes 
a fermionic direction. Altogether, the component fields in the matter multiplet are 

Z G ft°(E,C 4|A ^£) 

p G im°(Z,C 4lAf ®C 2 ®K% 2 ) (3.3) 
Y G 0°(E, C 4 '-^ (8) Ky, (g> £ _1 ) . 

Of course, the Z represent the pullbacks, now to E, of homogeneous coordinates of CP 3 |A/ ". We 
will sometimes decompose Z 1 as Z 1 = (Z^\x A ) = (p a ,X a \x A ) i n t° its bosonic and fermionic 
components. 
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3.1.1 The infinity twistor 

In order to write an action for these fields, we must pick some extra data. This is a choice 
of constant, graded skew symmetric bi-twistor Iu, known as the 'infinity twistor'. Graded 
skew-symmetry means that T u = -(-l)\ IJ \lj U where \IJ\ = 1 if both I and J denote 
fermionic directions, and zero otherwise. For any two twistors Z\ and Z2, we will usually 
denote TijZ\ Z2 by (Z±, Z2). 

Projective twistor space carries a natural action of SL(4|AA; C) acting as linear transfor- 
mations on the homogeneous coordinates. This is the complexification of (the double cover 
of) the space-time AA-extended superconformal group. The role of the infinity twistor is to 
break conformal invariance and determine a preferred metric on space-time [62-64]. Specif- 
ically, if X ah = Z^Z^ are homogeneous coordinates for the bosonic part of the twistor line 
Z\Zii representing a point x in space-time, then 

, 2 e a bcddX ab dX cd 

is the space-time metric. According to this metric, lines in twistor space that obey I-X = lie 
'at infinity' in space-time. The fermionic-fermionic components Tab were examined in [65, 66] 
where it was shown that they likewise define a metric on the R-symmetry group. A non-trivial 
Tab thus corresponds to gauging the R-symmetry. 
For definiteness, we will make the choice 

t . n \ 

(3.5) 



2"ab 








Tab 



Tij-- 

where the even-even components Z a t, and odd-odd components Tab are given by 7 

T ah = e °J and Tab = ^^ab, (3.6) 

respectively. This infinity twistor is non-degenerate. Its inverse is T IJ /A, where 



(3.7) 



T IJ defines a holomorphic Poisson structure T IJ A of homogeneity —2 on twistor 
space. It will play an important role in the vertex operators. 

In (3.6)-(3.7), A is a constant of dimensions (mass) 2 . The powers of A can be understood 
as follows. Since x has dimensions (mass) -1 and A has dimensions (mass) a , the incidence rela- 
tions fi a = x aa X a show that [i has dimensions (mass) - 2. Similarly, the space-time fermionic 











T IJ = 


Ae Q/3 







V 


VA5 AB 



7 The fact that Xu with lower indices contains e a/3 with raised indices originates with our conventions that 
the bosonic twistor components Z a are written as A Q ). 
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coordinate i? has dimensions (mass) 2 ; so the twistor space fermionic directions x are dimen- 
sionless. The powers of A ensure that both 

XjjZ'dZ 3 = Ae^d// + e a P\ a d\p + VA5 AB X A d X B (3.8) 

and the Poisson structure 

^ii d d A z d d , d d rr c-ar d d , 

x w A w - ' °V A 9? + A "'S A s\, + W (3 ' 9) 

have homogeneous dimension (mass) . Bosonically at least, this dimension is important 
in ensuring that (3.4) indeed has dimensions (mass) -2 as expected for a space-time metric. 
Recall that the n-particle g-loop gravitational scattering amplitude comes with a factor of 
K 2g-2+n ^hese dimensions must be balanced by a total of 2g — 2 + n powers of X. 

With the choice (3.6), the incidence relations show that the space-time metric (3.4) 
becomes 

2 rj^dx^dx" 

ds = -jt . 9N9 (3.10) 

(1 + Ax 2 ) 2 

where r]^ u is the flat metric. This is the metric of (complexified) AdS4 with cosmologi- 
cal constant A. Similarly, with Tab = V^$ab the SL(A^;C) R-symmetry is broken to 
SO(A/";C). Thus, with the choice (3.6) of infinity twistor, our model will describe (subject 
to an appropriate reality condition) SO (A/") gauged supergravity on an AdS4 background. It 
is straightforward to introduce an arbitrary gauge coupling for the gauged R-symmetry by 
rescaling Tab — > gv'A^AS for some dimensionless coupling g. In section (4.1) we shall take 
the limit A — > (with g remaining fixed) so as to compute the flat space S-matrix of ungauged 
supergravity. Until then, we set A = 1 and g = 1 so as to lighten the notation. 

3.1.2 The action 

Having chosen our infinity twistor, the action for Z is simple to state. We have 8 



Si = ^ / d l \ 2 z(Z,dZ) 



(3.11) 



where p aI = e at 'p^. Notice that the charge +2 of (Z, dZ) balances the charge —2 of d 1 ' 2 ^. In 
writing this action, we let d denote the (covariant) Dolbeault operator acting on sections of 
the appropriate bundles — since the (0,2)-part F 0,2 of the curvature of any bundle vanishes 
trivially on restriction to a Riemann surface, we can always work in a gauge in which the 
(0,l)-form part of any connection is the usual <9-operator. (This applies equally to the cs 
manifold X.) 



3 In computing the component expansion of this action, we use the conventions J d 2 8 8 a 8 a = 2 where 



9 a = e ab 8 a 8 b . Note also that 6 a 6 b = -\e ab ( 
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Clearly, to have sensible kinetic terms for all components of Z, it is important that the 
infinity twistor ( , ) be totally non-degenerate, both in the bosonic directions and fermionic 
directions. This motivates our choice (3.6). When we come to take the flat space limit in 
section 4 so as to describe scattering amplitudes, the infinity twistor necessarily becomes 
degenerate. We shall then need to rescale the fields so as to remove the dependence on the 
cosmological constant from the action, at the cost of including it in the definition of the 
Z(z, 6) supermultiplet. 

The KZ-system 9 

S Y z = ^j^Y I dZ I (3.12) 

was a key ingredient of Berkovits' twistor string [56]. Here, as there, performing the path 
integral over Y will lead to the constraint that Z be a holomorphic section of C 4 '-^ ® C. 
Thus, on-shell, Z describes a holomorphic map to PT. Berkovits' twistor string describes 
non-minimal N = 4 conformal supergravity [67-69], at least at tree- level, as does Witten's 
original model [12]. The M = 4 (conformal) gravity multiplet is not self-conjugate under CPT 
transformations and to build a CPT invariant theory we must use two separate supermultiplets 
that are exchanged under CPT. These two multiplets have a very different character on 
twistor space. The multiplet containing the positive helicity graviton is described locally by 
a vector field V, whereas the multiplet that contains the negative helicity graviton is instead 
described locally by a 1-form B 10 . Correspondingly, in Berkovits' twistor string the two 
conformal gravity multiplets are represented by the worldsheet vertex operators V I (Z)Yj 
and Bj(Z) dZ 1 , respectively [68]. 

Conformal gravity, being a fourth order theory, contains twice as many on-shell degrees 
of freedom as Einstein gravity. If we wish to extract the Einstein supergravity multiplets 
from the vertex operators of the Berkovits twistor string, we should require that 

= (l IJ djh) y 7 and B = (Z, dZ) (3.13) 

for some (local) functions h(Z) and (f>(Z) of homogeneities +2 and —2, respectively. (See [24, 
62, 64, 66] — or section 3.4 below — for further details.) One of the challenges to be overcome 
in constructing a twistor string for J\f = 8 supergravity is to understand how to unify these 
'vector field' and 'one form' vertex operators as part of a single CPT self-conjugate M = 8 
multiplet. 

Although it is premature to discuss the spectrum of our model at this point, (3.2) already 
contains a small hint of the solution: the fields Z and Y are unified into a single worldsheet 

Here we have used the infinity twistor X to lower the index on Y, so that Yi = 1ijY j . With a non- 
degenerate infinity twistor, this is harmless, but in the flat space case this operation must be done with care. 
See section 4.1. 

10 More precisely [68], in JV = 4 conformal supergravity V represents an element of H 1 (¥T, Trt) and defines an 
infinitesimal deformation of the complex structure of a patch of twistor space that preserves the holomorphic 
section D 3|4 Z of Ber(PT). The conjugate field B plays a role similar to that of the heterotic £?-field. Its 
curvature H = dB represents an element of _ff 1 (PT, 
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supermultiplet. Thus, from the perspective of X, there is no fundamental difference between 
the two types of vertex operator in (3.13). 

3.2 BRST transformations 

To complete the specification of our theory, we need to choose a BRST operator. This will 
be based on the symmetries of X that act trivially on S, as discussed in section 2.1. 

Consider the following three sets of transformations of the component fields. Firstly, 

b v Z l = e a p{ , Sipi = EaY 1 , dyY 1 = (3.14) 

with fermionic parameters e a , secondly 

52Z 1 = ^Z 1 , 5 2 p[ = -K a b p , 62V 1 = -^Y 1 (3.15) 

with bosonic parameters K b , and finally 

5 3 Z T = 0, 5 3 pi = -^Z 1 , ^y^l^pi (3.16) 

with fermionic parameters e a . These transformations represent the actions of a local symme- 
try of X on the matter multiplet Z. On a local coordinate patch U C X, we may take the 
parameters (e, k, e) to be constant. In this case the action Si is invariant under (3.14)-(3.16) 
when restricted to U. However, the spins and charges of the component fields mean that if we 
wish to make sense of these transformations globally over S, then e and e must transform non- 
trivially on overlaps. Specifically, we must have e a 6 77fi°(£, £<g)K E 1//2 ), K a h G f2°(£, O) and 
e a G 7IO°(S,£- 1 ® K^ 1/2 ), so that they fit together to form a supermultiplet in S1°(X,£>). 
In particular, to treat (3.14)-(3.16) globally over E, the parameters e and e must depend 
(smoothly) on the worldsheet coordinates. Thus these transformations must necessarily be 
gauged. 

3.2.1 The ghost multiplets 

With non-constant parameters, the matter action is not invariant under (3.14)-(3.16). To 
remedy this, and to treat the transformations as redundancies, we follow the usual procedure 
of introducing ghosts. Since the above transformations reflect the actions of S-preserving 
symmetries of X as studied in section 2.1, we introduce a ghost multiplet 

Cenn°(X,V) (3.17) 

in the parity reverse of the parameter multiplet. We declare C to have ghost number n g h = +1. 
As with the matter field, we can expand C in terms of components as 

C a (z, 6) = 7 a (z) + 9 b N a b (z) + -9% u a {z) , (3.18) 

where again a = 1,2. Recalling from section 2 that UV ^ C 2 <g> 1/2 ® C, we see that the 
component fields 7° are a pair of are (— |,0)-forms of charge +1. They are bosonic ghosts 
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for the fermionic parameters e a in (3.14). The bosonic fields v a are likewise 0)-forms of 
charge —1 and are ghosts for fermionic parameters e a in (3.16). Finally, N b are four fermionic 
ghosts corresponding to the GL(2; C) transformations of the rank 2 bundle C 2 (g) C They 
are scalars on £ of charge 0. We shall often find it convenient to separate this GL(2; C) as 
GL(1;C) x SL(2;C). Accordingly, as in (2.17) we write 

W b = l -5 a b n + n\ (3.19) 

where n a b is symmetric and traceless and n = tr(N). Note that since n is the (anticommut- 
ing) gauge parameter for the determinant of C 2 <8> £, the appropriate parameter for gauge 
transformations of C itself is n/2. This explains various factors of \ that appear in the BRST 
transformations below. 

We also introduce an antighost multiplet 

B € 77ft°(X, Ber(X) ® £> v ) (3.20) 

of Tigh = — 1 that is conjugate to C. This may be expanded as 

B a (z, d) = na{z) + b M ab (z) + ^9 b 6 b p a {z) . (3.21) 
bmce i7Ber(X) ® £> v S C 2 ® i-sT+ 1/2 ® £ > 

we see that the two bosonic fields /x a are (^, 0)-forms 
of charge +1, the fermionic antighosts M a b are uncharged (l,0)-forms, and finally f3 a are a 
pair of bosonic (|,0)-forms of charge —1. As with the ghost N, we shall often separate the 
antighost M into its GL(1;C) and SL(2;C) parts, writing 



M ab = e ab m + m ab (3.22) 



with m ab symmetric and traceless. 
The ghost action is simply 



S 2 = i I A l \ 2 zB a dC a 
27T Jx 

= ^ ! Ml a + m ab Bn ab + mdn + ^ a du a 
^ Jt, 



(3.23) 



Except for their non-trivial charges under C, the /?7-system is just two copies of the usual 
/?7 system of the RNS superstring, the MN-system is the standard system associated to 
fixing the GL(2; C) symmetry of V, and the /zzz-systems corresponds to fixing supergauge 
transformations associated to gauginos. The non-trivial charges of these ghost fields mean 
that the gravitinos and gauginos that they fix are also charged. 

The above behaviour is perhaps reminiscent of a GL(2; C) gauged supergravity on the 
worldsheet. However, because we are only gauging those symmetries of X that act trivially 
on the bosonic Riemann surface S, we are not actually considering worldsheet gravity itself. 
Correspondingly, our ghosts live only in the subgroup 770° (X, P) of i7f2°(X, TX) and there 
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is no (fundamental) fermionic foc-ghost system. It may seem strange to have gravitinos (albeit 
non-propagating ones supplanted by the /37-system) but no graviton. Usually in supersym- 
metry, this is not allowed because the structure of the supersymmetry algebra {Q, Q^} = P 
forces us to gauge Poincare transformations if we gauge the supersymmetry. In the present 
case it is possible to have gravitinos without gravitons ultimately because the distribution T> 
is integrable and {T>, T>} 7$ T£. 

3.2.2 The BRST operator 

In the presence of the ghosts, the transformations (3.14)-(3.16) are replaced by BRST trans- 
formations generated by the operator 

Q = \j d 112 ^ (C a (Z, D a Z) + B a {C, CT ) , (3.24) 

where the derivative D a = d/d9 a and where {C, C} a = 2C b D b C a denotes the anticommutator 
- i.e., the graded Lie bracket on the distribution T>. The integral is to be taken over a real 
1-dimensional cycle T C S as well as over the fermonic directions. Clearly, Q is a fermionic 
operator of n g h = +1, and the spins and charges of the fields and measure d 1 ' 2 ^ combine to 
ensure that Q is a scalar of charge zero under GL(2;C). It is also important to notice that 
the BRST operator depends on our choice of infinity twistor ( , ). 

Performing the integrals over the anticommuting coordinates 6 a , (3.24) may equivalently 
be written as 11 

Q = (fr(Y,Pa) + ^ a (Z, Pa ) + - (Y,Z) - ~n a V,p 6 ) 

J 1 1 1 v (3.25) 

U a _L n a .-v b \ 4- (--u a + n a ,u b ) 4-rrw a i,_ -m. L ( n( a n fe ) c 4- ^/Ky 6 )" 



+ /3 a (- 7 a + n a 67 °) + p a {--v a + n>°) + mrfv a - m ab [n^ c n°> c + ^ a v 

in terms of the component fields (3.2), (3.18) & (3.21). When acting on the matter multiplet 
2(z,6), this operator generates the transformations 

sz 1 = Vtf + ¥ 

Sp'^jaY 1 + ^aZ J -n a b pi (3.26) 

5Y 1 = 1 v a p I - U Y 1 
2 Ha 2 

generalizing (3.14)-(3.16). Similarly, the BRST transformations act as 
Sj a = *j 7 a + n a b j b bv a = -% a + nV" 

<5n = v a la 5n ab = n( a c n h)c + ^ a v b) 



(3.27) 



xl Our conventions are that /' a <? 6 ' is the symmetrized product \{.f a g b + f b g a )- Two-component indices 
a, b, are raised and lowered using the SL(2)-invariant anitsymmetric tensor e a b- In particular, for bosonic 
fields such as 7 and v, ^ a v a = eabj"^ = — ebai jh ^ a = — v h ^b- Pairs of anticommuting fields would have an 
extra minus sign. Some care has been taken to ensure the numerical factors in (3.25)-(3.28) are correct and 
compatible with the numerical factors in the matter and ghost action. 
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on the ghost multiplet. These transformations directly reflect the structure of the alge- 
bra (2.15)-(2.16). For example, we see that 7° transform in the fundamental of SL(2; C) and 
have charge +1 under C, while v a is again in the fundamental of SL(2; C) but has charge 
— 1. The unusual factors of ^7 in the transformation of n and n ab come from the fact that 
{vi,V2} = R'12 m (2.15)-(2.16). Finally, the BRST transformations act on the antighost 
multiplet as 

1 n 

5m=l((Z,Y)-/3 a7 a + f i a v a ) 

1 (3.28) 
5m ab = -(p a , p b ) - 2n (a c m fe)c + fi {a j b) + p {a u b) 

S/3 a = (Pa, Y) - ~n/3 a - n b (3 b - mu a + m ab u b , 

giving the currents conjugate to each symmetry. The transformations (3.26)-(3.28) have been 
checked to be nilpotent and to be symmetries of the full action Si + S2. Of course, the 
statement that Q 2 = is subject to potential anomalies — so too is the very definition of the 
composite Q operator itself. We shall investigate these anomalies in section 3.3 below. 

BRST invariant configurations may be found by setting the fermionic fields to zero and 
asking that they remain zero under a BRST transformation. In this regard, the most dan- 
gerous looking transformation is 

5pi = KaZ 1 (3.29) 

which potentially forces Z to vanish, ruining the interpretation of our model as a map to 
projective twistor space. Even if only some components of the supertwistor Z 1 were forced 
to vanish, this would still place intolerable restrictions on the map and destroy any chance of 
the model describing (non self-dual) gravity. Of course, the resolution is that in fact v = 0. 
To see that this is so, notice that since deg£ > 0, the field v a G 0°(S, K^/ 2 (g> C^ 1 ) has no 
zero- modes (at least generically, and always at g = 0). Thus the path integral 

Vp ex V Q- J^p a du a ^j (3.30) 

over the non-zero modes of the conjugate fields p imposes the constraint v = 0, render- 
ing (3.29) harmless. Likewise, the transformation 5m ~ (Z, Y) + ■ ■ ■ is tame because Y 
vanishes on-shell. 

To summarize, our model contains just three field multiplets 

z e n°(x,c 4| ^®£) 

c e nn°(x,v) (3.31) 

B € 77fi°(X,Ber(X) ®£> v ) , 
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and is defined by the action 

S = ^f x dM2z {\{Z,dZ) + B a dcA (3.32) 

and the BRST charge Q of (3.24). This is enough to describe perturbative Af = 8 super- 
gravity, also allowing for gaugings and conformally flat backgrounds such as AdS4. The only 
structures involved are the choice of infinity twistor ( , ) and the structure of X as a complex 
supermanifold. As promised in section 2.1, zero modes of C represent (parity reversed) global 
automorphisms of X that act trivially on E. Similarly, as in section 2.2, parity reversed zero 
modes of B are Serre dual to i? 1 (X, V), the tangent space to the moduli space of X as a 
bundle over E. 

3.3 Worldsheet anomaly cancellation 

The worldsheet theory we have defined is chiral — the matter and ghost kinetic terms each 
involve only the worldsheet Dolbeault d operator — so it is potentially rife with anomalies. 
We now investigate these, showing that all (local) anomalies cancel when Af = 8. 

We first compute the anomalies in the worldsheet gauge theory. Consider first the 
GL(1;C) transformations associated to the non-trivial line bundle C. On the two dimen- 
sional worldsheet, this anomaly is governed by a bubble diagram with the charged chiral 
fields running around the loop. It is thus determined by the sums of the squares of the 
charges of the fields, weighted by a sign for fermions. The GL(1) charged fields are the YZ 
system, giving a contribution a = (4 — Af) to the gauge anomaly, the two /?7 systems each 
giving a = 1, and the two \iv systems that also contribute a = 1 each. Altogether we have 

QGL(i) = (4-A0 + 2 + 2 = 8-AA (3.33) 

so the GL(1; C) gauge anomaly vanishes if and only if Af = 8. 

If the SL(2;C) bundle has non-trivial second Chern class, there is a further potential 
gauge anomaly. From the matter fields, only the pp-system transforms non-trivially under 
SL(2 C), in the antifundamental representation. The bosonic /J7- and //^-ghosts transform in 
the fundamental (or antifundamental), while the fermionic m a ;,n a6 -system is in the adjoint. 
The anomaly coefficient is thus 

a S L(2) = -^(4-A0tr F (A fc ) + 2tr F (t k t k ) -tr adj (A fc ) , (3.34) 

where t k denote the generators of SL(2;C), in the representation indicated by the subscript 
on the trace. (A sum over k is implied.) Note that the pp contribution has a symmetry factor 
2 since they are their own antiparticles. Writing ti^(t k t k ) = C*2(R) dim(R) in terms of the 
quadratic Casimir of the representation, (3.34) becomes 

a S L(2) = y x 2C 2 (F) - 3C 2 (adj) = ^(M - 8) , (3.35) 

12 Recall from (2.12) that Ber(X) = £ 2 and that d 1 ' 2 z is a holomorphic section of the trivial bundle Ber(X) <g> 
£~ 2 O. 
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where we used the SL(2;C) quadratic Casimirs 62(F) = | and C2(adj) = 2. The worldsheet 
GL(2; C) gauge theory is thus completely free from local anomalies when J\f = 8. 

In addition, we can compute the total Virasoro central charge 13 . From the matter fields, 
the YZ-system contributes central charge c = 2(4— AA), twice the (complex) bosonic dimension 
of the non-projective target space, minus twice the fermionic dimension. The pp-system are 
spin 2 fields on the worldsheet of opposite statistics to Z. As in the RNS superstring, they 
contribute a further (4 — M) to the central charge. From the ghosts we have two bosonic 
/37-systems each contributing c = +11 as in the RNS string, while fermionic ghosts for the 
gauge system contribute c = —2 x dim(GL(2)). Finally, because the spin | /i^-systems are 
bosonic, they contribute c = — 1 each. The total central charge is thus 

c = 3(4-Af) + 22-8-2 = 3(8-A0 (3.36) 

and vanishes if and only if Af = 8. There is in addition a potential mixed GL(1;C)- 
gravitational worldsheet anomaly b which equals 8 — N and again vanishes when J\f = 8. 
Since tr R (t fc ) = for SL(2; C), there is never any mixed SL(2)-gravitational anomaly. 

For an alternative (though equivalent) view of things, performing the path integral over 
the non-zero modes of all fields leads to determinants of (^-operators. These 9-operators 
act on sections of bundles as appropriate for the charges and spins of the fields, and the 
determinants appear in the numerator or denominator according to whether the fields are 
fermionic or bosonic. The resulting chiral determinants are not functions, but form a section 
of a determinant line bundle over the moduli space of the gauge theory (and, in principle, 
the complex structure of the worldsheet). In order to make sense of the determinants as 
these moduli vary, we must find a flat connection on this determinant line bundle. A natural 
connection was provided by Quillen [70]. Its curvature F is computed by the Freed-Bismut 
formula [71, 72] 

T = / Td(TS) A Ch(£) , (3.37) 

where the bundle E depends on the fields in question. Letting x denote the first Chern class 
of TS, y denote ci(£), and Q denote the SL(2;C) bundle, we find 

t = ii Jx x Ax + ^ L x Ay + ^ L y Ay ~ s L c2(g) (3 ' 38) 

where cigl(i)! a SL(2)> D ! and c are the anomaly coefficients computed above. 

These potential anomalies would also afflict the BRST charge Q in (3.25), since it is a 
composite operator. For example, the terms proportional to the ghost field n each contain 
potential short distance worldsheet singularities. However, the coefficient of this short dis- 
tance singularity in the combination n ((Y, Z) — f3 a ^ a + [i a v a ) is the same anomaly coefficient 



13 Since our theory does not involve worldsheet gravity — although there is a worldsheet gravitino associated 
to the /37-system — the role of this anomaly is not completely clear to me. Its vanishing nonetheless seems 
significant. 
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a = (4 — J\f) + 2 + 2 as before, so Q is well-defined when M = 8. Similarly, the terms propor- 
tional to n a f) are sensitive to any anomaly in the SL(2; C) gauge symmetry. Finally, there are 
short distance singularities that potentially obstruct Q 2 = 0, so that Q could not be used as 
a BRST operator. It is left as an exercise to show that these again cancel when J\f = 8. 

3.3.1 Zero modes 

The absence of gauge anomalies means the path integral over the complete set of non-zero 
modes of all fields is well-defined, providing a section of a determinant line bundle over the 
moduli space whose Quillen connection is flat. We now examine the properties of the zero 
modes of the fields. 

Consider first the charged fields in our theory. For a generic chiral /?7-system, with 7 
taking values in some vector bundle E, Serre duality and the Riemann-Roch theorem give 

n 7o- n /?o = h°(Z,E)-h 1 (Z,E)= I Cl (det#) + Jci(TS) . (3.39) 
In the case at hand, writing d for the degree of C, we find 

nzo " n Yo = 4|AA x (d + 1 - g) (3.40) 
for the bosonic and fermionic components of the yZ-system at genus g, 

n 70 - uft, = 2(d + 2 - 2g) (3.41) 
in total for the two /37-systems in our theory, and a total of 

TVo " TVo = 2 d (3-42) 

for the two /^-systems 14 . 

The first important consequence of this calculation is that the path integral measure 
V(Z, Y, (3, 7, fi, v)q over these zero modes has net charge 

(4 - AT)(d + 1 - g) + 2(d + 2 - 2g) + 2d = (8 - Af)(d + 1 - g) , (3.43) 

where we recall that by Berezin integration, the integral form d# scales oppositely to 6 for a 
fermion. Thus, when J\f = 8, the zero mode path integral measure provides a top holomorphic 
form on the moduli space of the theory. 

However, while the total charge of all the zero modes cancels, there are selection rules 
associated to the zero modes of the individual fields. Generically, when d is sufficiently larger 
than g (and always at g = 0) the Kodaira vanishing theorem asserts that the above indices 
are entirely due to the positively charged fields. Let us examine the consequences of these 
selection rules, concentrating on this generic case. 

14 Here it is the antighost /1 that has zero modes, at least generically. 
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To begin with, we have (d + 1 — g) zero modes of Z A = \ A for each A = 1, . . . , Af 
running over the fermionic directions of twistor space. These fermionic zero modes cause the 
path integral to vanish unless they are saturated by insertions from the vertex operators. 
Exactly as in the original twistor string [12, 56], this leads to a relation between the degree of 
the curve Z(E) C CP 3 '"^ and the allowed helicity sector represented by the vertex operator 
insertions. We will find that the vertex operators describe an M = 8 gravity supermultiplet, 
in conventions where the positive helicity graviton is at order (%)° and the negative helicity 
graviton is at (x) 8 - Thus we find the usual relation 

d = k+l + g (3.44) 

between the MHV level k ('number of negative helicity gravitons, minus 2') and the degree 
of the curve and genus of the worldsheet 15 . 

Now we consider the zero modes of the bosonic fields. We have d + 2 — 2g zero modes for 
each of the two components of 7 a and d zero modes of each of the two fj, a s. Of course, we also 
have (d + 1 — g) zero modes for each bosonic twistor component. In the absence of insertions 
that depend on these zero modes, the integrals over 7 and \i would diverge, rendering the 
path integral ill-defined. As in the RNS string, we will find that the vertex operators have 
(^-function support in these fields, giving a meaningful path integral. For now, recall from the 
introduction that in the flat space limit, an n-particle g-loop amplitude with n± gravitons of 
each helicity is a monomial of degree u_ — 1 + g in the infinity twistor ( , ) as a form, and 
n + — 1 + g in the infinity twistor [ , ] as a Poisson structure 16 . With the aid of (3.44), we 
may rewrite these numbers as 

n_ — 1 — g = d 

n+-l-g = n-(d + 2-2g), 

respectively coinciding with the number of zero modes of each component of fj,, and n minus 
the number of zero modes of each 7 component. 

The remaining fields are the worldsheet spinors p, which generically have no zero modes, 
and the ghost system associated with gauging the GL(2; C) transformations. These ghosts 
certainly do have zero modes. However, at g = we will content ourselves to treat these by 
simply 'dividing by vol(GL(2))' - the path integral will lead yield a form that is invariant and 
basic with respect to a natural GL(2) action, and we descend to the moduli space. While this 
suffices to recover the g = scattering amplitudes of [32], it is really too naive. We discuss 
this further in section 5. 

lj The genus h of the image curve Z(E) obeys h. < g. In particular at MHV level for g = 1, the worldsheet 
must double cover a twistor line, branched over four points. 

16 This monomial behaviour of course refers to the amplitude when written in twistor space. It becomes 
obscure on momentum space because the transformation from twistors to momenta itself involves the infinity 
twistor. Recall also that I is degenerate in the flat space limit. 
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(3.45) 



3.4 Vertex operators 

In this section we construct the vertex operators representing BRST cohomology classes. 
These will correspond to a linearized N = 8 supergravity multiplet. We also construct 
picture changing operators required to fix zero modes of the bosonic antighosts. 

3.4.1 The M = 8 supergravity multiplet 

The odd supervector field V of (2.14) generates a global holomorphic automorphism of X 
when V 6 H (X, X>). In the generic case that d ^> g, only the lowest two components v and 
R in the superfield expansion of V can be globally holomorphic, with 

ve#°(£,£>) and R € #°(£, End(C 2 ® £)) . (3.46) 

The fermionic symmetries of X — > S correspond to zero modes of the ghosts 7, with each 7" 

— 1/2 

being one of the d + 2 — 2g holomorphic sections of Kj, (g) £, while the bosonic symme- 
tries correspond to the zero modes of the ghost N a b that are constant. To obtain a moduli 
space whose (virtual) dimension is non-negative, and hence to have a well-defined ghost path 
integral, we must fix these zero modes. 

In our case, the odd vector field v = v a d/d6 a s generates (smooth) translations of along 
the fibres of X — > S. To fix the associated odd automorphisms of X we pick points p% £ T, 
and demand that the translations act trivially at these points. In the path integral, these 
translations are represented by the ghosts j a , so we can force the translation to be trivial at 
some pi by inserting 17 

S\ 1 ) = 5( 1 1 )S( 1 2 ) (3.47) 

at this point. Each such constraint reduces the dimension of the space of automorphisms 
by 1, so generically we need to pick (at least) d + 2 — 2g points to remove all the global 
odd automorphisms. The resulting <5-functions absorb the 7 zero modes, rendering the path 
integral meaningful. 

In the usual case of the RNS superstring, the vertex operator would also include a factor 
of the fermionic c ghost instructing us to quotient the path integral only by those diffeomor- 
phisms of the bosonic Riemann surface X that act trivially at the pi. For us however, since 
we only quotienting by diffeomorphisms generated by sections of T> C TX, there are no c 
ghosts. If we do not wish our answer to depend on the choice of pi G S, we must integrate 
over them 18 . Because 7 a has spin —\ and charge +1 under C, the operator 5 2 {^) should be 
interpreted as a (l,0)-form of charge —2. So the simplest type of vertex operators are 

O h = I 5 2 { 1 )h(Z) , (3.48) 

17 The fact that our vertex operators can depend only on 7, and not on derivatives of 7, is determined by 
the requirement that the path integral measure over all fields does actually descend to a measure — or top 
holomorphic form — on the moduli space of X — > T>. See [55] for an explanation in the context of RNS 
superstrings. 

18 See however the discussion in section 5.1. 
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where h(Z) is a (O,l)-form on £ of charge +2. This integral is to be taken over £ at 6 a = 
0. These vertex operators are closely analogous to Neveu-Schwarz vertex operators in the 
superstring. 

The field h is a twistor representative of an M = 8 supergravity multiplet, pulled back 
to S. Writing Z 1 = (Z a |x^) for the bosonic and fermionic components of the twistor, we can 
expand h as 

h{Z\ X ) = h(Z) + X A MZ) + \x A X B aAB(Z) + ■■■ + {xfh{Z) (3-49) 

where the coefficient of {x) p is a (0, l)-form on twistor space of homogeneity 2 — p. Via the 
linearized Penrose transform [73-75], these states correspond to massless fields with one boson 
of helicity +2, 8 fermions of helicity +|, 28 gauge fields of helicity +1 and so on until we 
reach the field h that corresponds to a negative helicity graviton. More precisely, the Penrose 
transform asserts that an on-shell, linearized N = 8 supergravity multiplet corresponds to a 
cohomology class in ^(PT, 0(2)), of which h is a representative. 

While Oh is to be inserted on a fixed section of X — ¥ £ (taken to be the zero section 
9 a = 0) , we also have related vertex operators Oh that are integrated over the entire worldsheet 
supermanifold X. As usual, these are obtained simply by replacing the factor of 5 2 {^) in (3.48) 
by the integration measure d 1 ' 2 ^ on X. We have 



O h = [ d^ 2 zh(Z) 

f dh(Z) . _ 1 b , j 2 h(Z) (3 - 50) 
i s dZ 1 2 f aPb dZ r dZ J ' 

Note that if /, J correspond to bosonic twistor directions, then p 1 and p J anticommute, and 
the e ab symbol ensures that the second term here is symmetric in /, J. Conversely, if /, J 
are fermionic directions, the expression is naturally antisymmetric in J, J. In particular, this 
means that (3.50) is well-defined as a composite operator, with no short distance singularities 
— the potential singularity in either the pp-system or the YZ-system are each proportional to 
the (graded) antisymmetric infinity twistor I, so the resulting derivatives on h would vanish. 

The first term in Oh also has a natural meaning in twistor space: since h represents an el- 
ement of ^(PT, 0(2)), and in canonical quantization of the action (3.11) we have Y ~ d/dZ, 
the term -j^jY 1 represents an element of i? 1 (PT, Tpx). It therefore describes an infinitesimal 
deformation of the complex structure of twistor space. The non-linear Penrose transform [62] 
asserts that performing a finite deformation of the complex structure of twistor space corre- 
sponds to turning on self-dual Weyl curvature in space-time. The holomorphic geometry of 
twistor space determines the conformal structure of space-time, so an arbitrary deformation 
of this complex structure, generated by an arbitrary V G H l (YY, 2ft), corresponds to an 
arbitrary self-dual solution of the Bach equations of conformal gravity. However, unlike in 
the original twistor string, the vertex operators we have found here are associated not to 



- 24 - 



arbitrary vector fields, but rather to vector fields 19 V£ = that are Hamiltonian 20 

with respect to the Poisson structure defined by X. Deforming the complex structure by a 
such Hamiltonian vector field ensures that the holomorphic Poisson structure is preserved, 
and hence the corresponding deformed space-time still has a preferred metric. This metric is 
a self-dual solution of the vacuum Einstein equations [62-64]. Extending this to the N = 8 
multiplet gives a BPS solution to the field equations of supergravity. 

As usual, the integrated vertex operator (3.50) can be added to the original action 



describing strings propagating on a background twistor space with deformed complex struc- 
ture d — > d + Vh- In the twistor string framework, deformations that are not self-dual are 
described perturbatively in terms of higher degree maps (worldsheet instantons). 

To summarize, freezing all the 7 zero modes requires us (in the generic case) to include 
d + 2 — 2g vertex operators of the form Oh, involving the basic twistor wavefunction h. The 
remaining n— (d + 2 — 2g) states are represented by integrated vertex operators Oh that 
involve not h itself, but rather its associated Hamiltonian vector field Vh- In the introduction 
we saw that the flat space tree- level scattering amplitude involves precisely n — (d + 2 — 2g) 
powers of the infinity twistor [ , ] (i.e., the infinity twistor as a degenerate Poisson structure). 
In particular, at g = the tree amplitudes (1.5) are monomials of degree n — d — 2 in [ , ]. 
We now understand that this fact has its origin in the odd automorphisms of X as the total 
space of a fermionic bundle over a fixed S. 

3.4.2 Picture changing operators 

The worldsheet supermanifold X also has moduli, even for a fixed Riemann surface E, coming 
from the freedom to deform the distribution T> C TX. As in section 2.2, for fixed S the tangent 
space to this moduli space is i7 1 (X, T>). Our twistor string knows about these moduli via the 
zero modes of the antighost multiplet B, which live in the parity reversed Serre dual group 
nH°(X, Ber(X) <g> £> v ) ^ H°(X, C 2 ® k]! 2 <g> C). In the generic case with d > g, the only 
components of the B multiplet to have zero modes are \x and M. 

Consider first the odd moduli space, associated to zero modes of the bosonic \x antighost 
in i?°(S, C 2 ®K^/ 2 ®£). To integrate over the odd moduli space we follow the usual procedure 
of the RNS superstring and insert 2/i°(E, K^/ 2 ® C) 'picture changing operators' 



1 We have used T to put the vector field index in the natural place. This is harmless when T is non-degenerate, 
and in section 4.1 we shall see it happens automatically in the flat space limit. 
20 This is the reason we denote the gravity multiplet by h. 




(3.51) 



T = 2 5 2 (n)S a S l 



(3.52) 
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where 

■^o = 2 ^ Pa ^> + 2 n/ifl ~ na ^ b + ™ 7a + mab76 ( 3 - 53 ) 
is the supercurrent obtained by taking the BRST transformation of p, a as in (3.28). These 
insertions correspond to 5-function wavefunctions for the gauginos associated to the \iv- 
system, and fix the p, a zero modes. See e.g. section 3 of [55] for a clear discussion of picture 
changing operators and their relation to fixing parity odd moduli. 

Unlike the usual picture changing operators of the RNS string, the operator (3.52) involves 
two copies of these currents because V has rank two. Because xS(x) = we can neglect the 
terms proportional to the \i antighost in these supercurrents. Similarly, at g = when the 
antighosts M a b have no zero modes (and there are no N insertions with which to contract), 
we can neglect their contribution to S a . Then the picture changing operator simplifies to 
become 

T = -S 2 (p)(Z,p a )(Z,p a ) (3.54) 

When g = we need (a minimum of) d such T insertions. 

The previous expression may be thought to be somewhat formal, because composite 
operator (Z, pi)(Z, p%) has a potential short distance singularity from the pp contraction. To 
regularize this, we point split the two (Z, p) factors and take a limit as they come together. 
In terms of a local coordinate z on U C S, one finds 

lim ((Z, Pl (z')) (Z,p 2 (z))) = lim x (Z(z'),Z(z)) , (3.55) 



where the factors of vdz' and \/dz arise since p is a spinor on S. Expanding the holomorphic 
field Z as Z^z') = Z 1 (z) + (z' — z) dzZ 1 (z) + • • • , we see that the pole from the pp propa- 
gator is cancelled by a zero from the antisymmetric infinity twistor, leaving us with a finite 
contribution 

lim ((Z, Pl (z')) (Z,p 2 (z))) =-{Z,dZ)(z). (3.56) 

z'— >z VP 

We now define T more precisely as the normal ordered operator 

T = S\p) Q :(Z,p a )(Z,p a ): - (Z,dZ)\ (3.57) 

in which the local contribution of the pp-system is explicitly accounted for. The normal 
ordering prescription : : is understood to mean that we do not consider contractions between 
the enclosed fields. 

Actually, since the potential short distance singularity cancelled in (3.55), we are free 
to think of T as in (3.54) without normal ordering. We must then remember to include the 
local contribution (3.56) when computing correlation functions involving these operators. In 
practice, this approach turns out to be somewhat simpler. 

Of course, we could have chosen to represent all the external states by the vertex operators 
Oh of (3.48), rather than use any integrated ones Oh- In this description, we would quotient 
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the path integral only by translations of S inside X that act trivially at n > d + 2 — 2g points. 
Since we are quotienting by fewer fermionic symmetries, the odd dimension of the moduli 
space increases and we must integrate over this larger odd moduli space. In the language of 
ghosts, the additional insertions of Oh provide extra factors of <5 2 (7). To compensate for these 
we should also construct picture changing operators for the /37-system, inserting n— (d+2 — 2g) 
of them so as to provide a top form on the odd moduli space. This is expected to be the 
correct approach if one wishes to obtain a detailed understanding of the compactification of 
this moduli space [55, 59]. It would be interesting to investigate this further. 

Finally, when g > 1 we also have zero modes of the fermionic antighosts M-ab- Insertions 
of these amount to constructing a top holomorphic form on the (bosonic) moduli space of 
holomorphic GL(2; C) bundles on E. For the (generic) case that this bundle is stable and 
g > 2, this moduli space has dimension 3(g — 1) + g and has been extensively studied [76-83]. 
We discuss it further in section 5.1. 

4 Scattering amplitudes in the flat space limit 

Our prescription for computing n-point worldsheet correlation functions in the g = twistor 
string is to consider the path integral 

n d ^ii T )=Cfi / ^)m^) n / dii2 ^(^) u r ) ■ ^ 

\j=l fc=d+3 1=1 I \j=l jT ' fc=d+3 X 1=1 I 

In this section, we will use this prescription to recover the flat space tree-level S-matrix of 
N = 8 supergravity in the form obtained in [32]. To do so, we will need to be able to handle 
correlation functions of /37-systems involving operators such as <5 2 (7). A clear explanation of 
how to achieve this was recently provided in [55] (see especially section 10). For convenience, 
the relevant points are summarized in appendix A. 

4.1 A degenerate infinity twistor 

To compute scattering amplitudes, we must take the limit as the cosmological constant A 
tends to zero. In this limit, the rank of the infinity twistor 
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we have been working with so far becomes non-maximal. In particular, in the flat space limit 
we must carefully distinguish between the infinity twistor as a form and the infinity twistor 
in its role clS ct bivector, since 

X LJ T JK = Adf ^0 (4.3) 

and so they are not equivalent. 



If we were to take the flat space limit naively, the matter action (3.11) would also become 
degenerate, with the kinetic terms for the \x a and \ A components of the supertwistor Z 1 
vanishing. To avoid this, before taking the limit, we relabel the fields as 
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(4.4) 



and include an overall factor of 1/A in the normalization of (3.11). In terms of the rescaled 
fields, the matter action becomes 

1 



Si 



2vr 



YtfZ 1 + pidp 1 



(4.5) 



and is independent of the cosmological constant. The ghost fields are unchanged. Having 
ensured the action remains non-degenerate, we can now freely take A — > 0, setting 21 
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We follow the standard convention that ( , ) denotes contraction by Ijj with downstairs 
indices, involving only the A part of Z, whereas [ , ] denotes contraction by X IJ with up- 
stairs indices and involves only the derivatives d/dp tangent to twistor space (or the As on 
momentum space). 

While the rescaled action is independent of the infinity twistor, the same cannot be said 
for the BRST operator 

(4.7) 



Q^t C r = \jd 1 \ 2 z{Z : D a Z) . 



As with the action, we first apply the recaling (4.4) with a non-degenerate infinity twistor 
and then take the limit A — > 0. The matter BRST charge becomes 



Qflat = f^Yip 1 + 7 2 [Y,p\ + \v\p,Z) + ^pjZ 1 

+ InYjZ 1 + \ (n 12 + n 21 ) p% + \n*{p,p) + \n^p,p\ 



(4.8) 



and the presence of the degenerate infinity twistor means that not all supertwistor components 
appear in all terms; for example, (p,Z) = p a X a while [Y,p\ = Yaf> a . A somewhat similar 
BRST operator occurs in Poisson sigma models, see e.g. [52, 53]. It would be interesting to 
explore this connection further. Again, the ghost BRST charge is unaltered. Similarly, in the 
flat space limit the vertex operators become 



5\ 7 )h(Z) 



\y dh ' 
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(4.9) 



L In the presence of an arbitrary gauge coupling g, we are taking the limit A — > 0, g\/A — > 0. See the 



discussion at the end of section 3.1.1. 
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for the external states and 

T = S 2 ( f i)(p,Z)p I Z I (4.10) 

for the picture changing operator. As promised, the integrated vertex operator Oh naturally 
depends on the Hamiltonian vector field [J^, ] associated to the infinity twistor as a Poisson 
structure. 

4.2 The worldsheet Hodges matrix 

We are now in position to recover the flat space gravitational S-matrix from the correlator (4.1) 
at g = 0. In this section we will show that the worldsheet Hodges matrix H in (1.6) arises 
from the correlation function of the matter vertex operators Oh and Oh- 

Firstly, we notice that the only insertions of 7 come from the 5-functions in the fixed 
section vertex operators Oh- These <5-functions serve to fix the integrals over the zero modes 
of 7 a , representing elements of H°(T>, Kj, 1 ^ 2 C). For each flavour 71 and 72, we expand 7 
as 

d+2 

7a(c) = r «i^*( <7 ) + non - zero modes , (4.11) 

i=l 

where the form a basis of the zero modes (written in terms of a homogeneous coordinate u— 
on the CP 1 worldsheet), and where the r a s are c-number constants. In [55], it was explained 
that for each flavour of 7°, the insertion of J- functions leads to 

\i=i Ifr 

See also the discussion in the appendix. Here, Y is the (d + 2) x (d + 2) matrix whose 
entries are Yy = Vn{<Jj). At genus zero, a basis of H°(T>, K^ 1 ^ 2 0(d)) is given by ^i(a) = 
(crd(r) _1//2 (T-i • • • cr- d + 1 . Computing this determinant and including both flavours, the path 
integral over the /37-system yields 

/d+2 \ d+2 

^fo))) = kl ... ad+2 |2 >< n M-;) , (4.i3) 

where \a% ■ ■ ■ 0-^+2! denotes the Vandermonde determinant 

ki - • • o-d.+2i = n (y) • ( 4M ) 

i<j 

ij6{l,...,d+2} 

This Vandermonde determinant is precisely the denominator factor (1.7) of det (H) in the 
introduction, here specialized to the case that we remove the first d + 2 rows and also the 
last d + 2 columns in computing a minor of the Hodges matrix H. (That is, we compute the 
(d + 3) rd principal minor.) 
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The (n — d — 2) x (n — d — 2) minor of H itself comes from the remaining part of the matter 
vertex operators. As a first step to understanding this, consider the n — d — 2 insertions of 
Oh and temporarily neglect the [Y, ^] terms. The remaining part of Oh is bilinear in the 
worldsheet spinors p and p. These are free fields on £. Since they have no zero modes, all p 
and p insertions must be absorbed by contracting them pairwise in all possible combinations. 
Insertions of p and p can be found both in Oh and in the picture changing operators T. 
However, with our degenerate infinity twistor, Oh involves only the a components of p, while 
T involves only the a components of p. The off-diagonal two point function (p a (a) p a (a')) 
vanishes, so the p a s from any Oh insertion can contract only with the p a s present in some 
other Oh insertion. Furthermore, the pieces Pa~§^~ + P A ^k in P 1 ' ~§^t in (4-9) may be ignored 
because there is always at least one unpaired p a left over in some Oh and at least one unpaired 
p a left over in some T, either of which causes the path integral to vanish. We conclude that 
we can consider the pp factors in Oh independently from those in T, and that we only need 
consider the a terms in Oh- 

This being so, consider the correlator 22 



C(a d+3 ,...,a n ) = / Y[ 
\fc=d+3 



_ d ( hdh k , y 



\fc=d+3 



d 2 h k 
dp A dpP 



P"Pc, >«—"n-" ' P a p & 

(4.15) 

coming purely from the pp-system. In terms of the homogeneous worldsheet coordinates a—, 
the two-point function of the p/9-system is 



(<7jd0i)2(<Tjd(7j 



(4.16) 



Using this propagator to perform all possible contractions in (4.15) yields 



C{cr a+3, • • • ,cr n ) 
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(4.17) 



where 
are 



j(0) 
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is the (d + 3) rd principal minor of the matrix HI^ ) whose elements 
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dZi ' dZj 
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dpi ' dpj 



for i ^ j 



(4.18) 



and zero on the diagonal. In (4.18) we understand that d/dZi acts on hi in the product (4.17), 
differentiating this wavefunction with respect to (the p a component of) its argument Z{<Ji). 
For example, if the external states are taken to be twistor representatives 



hi(Z((Ji)) 



^ 8 2 (Xi - Si\(ai)) exp (sip a (<Ji)\ ia + SiX A ((Ji)riiA 



(4.19) 



22 In this expression /i denotes the bosonic twistor component, not the antighost, as should be clear from the 
indices. 
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of momentum eigenstates, then — > Si-Sj[ij]/(ij). 

U(o) ig n ot quite the full worldsheet Hodges matrix (1.6), because the diagonal elements 
vanish. To recover (1.6) in its entirety, we now consider the additional effect of the [Y, ^1 
terms in the integrated vertex operators (4.9). With the degenerate infinity twistor, only the 
a component of Y is present here. This Y& cannot contract with any Z in T, either because 
there is no short distance singularity or because we would again be left with an unpaired p. 
However, the Yq, from any given Oh may contract with the Zs in the wavefunctions h(Z) in 
any of the remaining Oh operators, and also in the 'fixed' Oh vertex operators. 
Since Y € 0°(S, K s ® 0(-d)) and Z <G fi°(E, 0(d)), the YZ-propagator is 



(Yr(cTi) Z J {pj)) 



§ J (°"idcr. 



d+1 

n 



(arj) 



(4.20) 



where the a ar in the product on the right are arbitrary. This product ensures that both 
sides have the correct homogeneity. It arises because d f is ill-defined if / is a (0,l)-form of 
homogeneity d > —1; we are always free to modify d f — )■ d f + g where g is an arbitrary 
holomorphic section of 0(d), since this is annihilated by d. We can fix a choice of propagator 
by specifying d+1 points at which d f vanishes. This is the role of the product in (4.20). 
Of course, any meaningful expression — such as the Hodges matrix — is independent of the 
choice of these points (see [26, 32, 33]). 

Using this propagator, with just one [Y, ^] insertion we have 
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(4.21) 

where in the first step we integrated out the p fields using (4.17). The term in braces in the 
last line is one of the diagonal elements (the (d + 3) rd diagonal entry) that we were missing 
from the full Hodges matrix. 

We now prove inductively that the sum of Y and pp insertions in each integrated vertex 
operator Oh means that the total correlation function assembles itself into the complete 
worldsheet Hodges matrix H. To start, this is certainly true when there are only 2 integrated 
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vertex operators (corresponding to an n-point MHV amplitude). For 



n 



Y, 



dhj 
dZ 



d 2 hi 



11 ^ 



J'=3 



nv 



I22 + 





I21 



mi2 




(cridai)(cr 2 da 2 ) \\K 



=1 



un inii2 
I21 H22 



(a 1 dai)(a 2 da 2 ) JJ/ii 



=1 



(4.22) 

where the first equality follows because we must either take the pp term at both sites or at 
none 23 (If there is only one Oh insertion — which occurs only for the 3-pt MHV amplitude 
— we are forced to take the Y contribution as in (4.21); the p fields cannot contribute at all.) 

To perform the induction, assume that the worldsheet correlator correctly gives the de- 
terminant of the full Hodges matrix when there are m— 1 Oh operators, for some value of m. 
Then from (4.21) 



Y. 



dhi 



n<\ n 

j=2 k=m+l I 



hi x P(m-l)x(m-l)| 



Y[hi n( fJ i dcJ i 



(4.23) 



i=l 



1 (m-l)x(m-l)| 



where the sum in each of the diagonal elements of the (m — 1) x (m— 1) minor 
also runs over site 1 , since the Y insertions in the Ohj s leading to this matrix may additionally 
contract with site 1. The Hodges factors on the right hand side of this expression can be 
written as the determinant of an m x m symmetric matrix with Hy = for j / 1. On the 
other hand, instead choosing the pp term at site 1 gives 24 
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8ZfdZ l j=2 fc=ra+l 
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H 23 






Him 


H 2m 







[[^[[(ffjdffj- 

i=i j=i 



(4.24) 

where the first row and first column represent the possible choices of contraction for the 
additional p and p insertions. As before, the sum in the diagonal entries of the (m— 1) x (m— 1) 
Hodges matrix should be extended to run over site 1. It is now clear that (4.23) & (4.24) 
combine to give the determinant of the full Hodges matrix appropriate to m insertions of Oh 
and n — m insertions of Oh- 

23 To lighten the notation, we have temporarily reversed our convention and taken the first n — d — 2 vertex 
operators to be integrated and the last d + 2 to be fixed. This corresponds to computing the leading principal 
minor of H. 

24 The lines in the matrix in (4.24) are simply to distinguish contributions from the new insertions at site 1 
from the previous inductive step. H is not a supermatrix. 
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Combining this with the factor (4.13) from fixing the 7 zero modes shows that the vertex 
operators (4.9) contribute 



/d+2 n \ . n 

n ^ n ) = / det/ ( M ) nwzfa)) 

\i=l j=d+3 / J i=l 

to the twistor string path integral, in the specific case that we choose to remove the same 
d + 2 rows and columns in computing a minor of the full n x n Hodges matrix (here chosen 
to be rows and columns 1 through d + 2). In [26, 32] we are actually free to compute 
any (n — d — 2) x (n — d — 2) minor of H, provided we divide by the corresponding two 
Vandermonde determinants. We can arrive at these more general representations by also 
allowing 'intermediate' vertex operators that involve a single £(7) and an integral over the 
other 6. That is, if we wish to compute a minor of the Hodges matrix involving different rows 
and columns, we should allow 

j dM(7i(<7)) KZ)\ di=0 and j d9 1 S( l2 (a)) h(Z)\ 62=0 (4.26) 

as well as Oh and Oh- The rows and columns that we remove from the Hodges matrix 
correspond to the independent insertion points of £(71) and £(72)- More generally, it should 
be clear that since the amplitudes depend on HI only through det (H), there is actually a very 
large amount of freedom in the elements themselves. It would be interesting to know if these 
examples considered in [84] can be realized on the worldsheet. Of course, since det'(H) is 
invariant under arbitrary permutations of all n external states, the minimal case considered 
in (4.25) is sufficient to recover the amplitude. Indeed, permutation invariance is now seen to 
be a consequence of the usual fact that it does not matter which vertex operators we choose 
to be 'fixed' and which 'integrated'. Thus we have recovered one of the main ingredients in 
the formula (1.5) for the tree-level gravitational S-matrix. 

4.2.1 Self-dual M = 8 supergravity 

In the following section we will show that the remaining, conjugate Hodges matrix comes from 
the d insertions of picture changing operators T. However, in the special case that d = - 
corresponding to constant maps to twistor space — X has no odd moduli and no Ts need be 
inserted. This case is worth investigating separately. 

Since Z(a) = Z for constant maps, we obtain the 3-point MHV amplitude 

<5 2 ( 7 )MZ) f S 2 (j)h 2 (Z) f d^zh 3 (Z)) 
v Jt, Jx I £4 27) 

V^ZMn{Z) A{h 2 (Z), h 3 (Z)} , 

where we divided by vol(GL(2; C)) in lieu of fixing the zero associated to the worldsheet gauge 
theory. This absorbs the integration over the three vertex operators over S and ensures the 
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remaining integral is taken over the projective twistor space. As usual, the braces { , } denote 
the Poisson bracket associated to the infinity twistor X as a Poisson structure X IJ -J^j A . 
Note that the Poisson bracket itself has homogeneity —2, while each hi{Z) has homogeneity 
+2, so (4.27) is well-defined on the projective twistor space. 

This 3-point MHV amplitude is especially important because it is the vertex of the action 

S sd = y"D 3|8 ZA (hAdh+ ^hA{h,h} \ , (4.28) 

evaluated on on-shell states. S s d is the twistor action for self-dual J\f = 8 supergravity and was 
first obtained by Mason & Wolf in [66] . At the linearized level its equations of motion say that 
h represents an element of H^'^iFT, 0(2)), as we have used in our vertex operators. At the 
nonlinear level, the equations of motion assert that the almost complex structure determined 
by d + {h, } is integrable. Once again, the fact that we deform the complex structure only 
by Hamiltonian vector fields ensures that we have a solution of self-dual Einstein gravity, 
rather than self-dual conformal gravity [62, 63]. S s d is clearly analogous to the holomorphic 
Chern-Simons theory 

SsdYM = J D 3|4 Z A tr yA A dA + A [A, A]^j (4.29) 

that describes self-dual M = 4 super Yang-Mills in twistor space [12]. Just as (4.29) is the 
string field theory of the perturbative open B-model, we can interpret (4.28) as the string 
field theory of our twistor string, restricted to constant maps. 

4.3 The conjugate Hodges matrix 

When d > we must also account for the picture changing operators T = 5 2 ([i)(pi, Z)p2iZ I . 
These will supply the conjugate Hodges matrix H v . 

In the recipe (4.1) there are no insertions of /i or v except for those in the picture changing 
operators, which just suffice to absorb the \i zero-modes. Recall that at g = a \i zero-mode 
is an element of H°(E, K% 2 ® 0(d)). Thus, at g = we can expand each \i a as 

fi a ((r) = (ada) 1/2 M O01 ... ad l o--i • • • cr-d-i + n0 n-zero modes , (4.30) 

where the Ms are constants. Following the discussion of the appendix, we find from (A. 13) 
that the <5-functions associated to freezing the odd moduli produce a factor 

(n ^(mm) ) = 7 - A ^ — ( 4 - si ) 

\t\ l m ntiM^)ki--^dl 2 

The final factor is the Vandermonde determinant (1.9) appearing in the denominator of 
det'(IHI v ) in (1.5), specialized to the case that we compute the first d x d minor of H v . Notice 
that when d = 1 the unique fx zero mode is Li(a) = (ada) 1 ^ 2 for each of the two antighosts ji\ 
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and Therefore, in the case of the MHV tree, the Vandermonde factor in (4.31) is replaced 
by unity. This was the prescription taken in [32]. 

The numerator of the conjugate Hodges matrix comes from the associated supercurrents. 
These are (p, A) piZ 1 = (p, A) (pap a + p a X a + PAX A ) ■ Again, the pp-system has no zero modes, 
so we must absorb all these insertions by contractions. Since (p, A) cannot contract with the 
a or A components of p, the only term in the bracket that can contribute is p a X a 25 . Using 
the propagator 

{p a {ai)f/'(a m )) = dcf 77—, 4 -32 

(Lm) 

as in (4.16), performing all possible p-p contractions yields 

( II to A > p ax ^)) = |H^d| x II (<n&n) , (4.33) 

\l=l I 1=1 

where EI^ xd | is the first d x d minor of the matrix with elements 



(A(.,)A(. m) ) for ^ m and H v = _ W^dA^» 
(lm) [cridai) 



To obtain this result, recall that contractions of worldsheet fermions lead to a determinant 
of a matrix whose Z-m th entry corresponds to a propagator from site I to m. The diagonal 
elements arise as in (3.56) since we are using the form of picture changing operator without 
normal ordering, so must allow contractions between p and p at the same site — see the 
discussion in section 3.4.2. Indeed, when d = 1 this is the only contribution. 

The off-diagonal elements of H v in (4.34) are exactly the same as those in the conjugate 
Hodges matrix (1.8). However, the diagonal elements in (1.8) and (4.34) appear to be different. 
Let us now show that (1.8) can be simplified so that it takes the form (4.34) 26 . To begin with, 
notice that (\(ai), A(<r m )) / (lm) is everywhere finite, since the potential pole is cancelled by 
a zero in the numerator. Now consider the diagonal term 27 

hi w ri toO^(M 

of (1.8), where we have included a factor of (aidai). The only possible poles in 07 come 
from the factors (a r l) involving the reference points. However, a key point in [32, 33] was 
that (4.35) was completely independent of these points (see [33] for a contour integral proof 
of this). Therefore (4.35) is actually holomorphic in 07. Furthermore, (4.35) is a scalar of 

25 Recall from section 4.2 that, with the degenerate infinity twistor of flat space, there could be no contribution 
from cross-contractions of ps in T with any p in Oh- 

26 I am greatly indebted to Lionel Mason for pointing this fact out to me, using a slightly different argument 
to the one given here. 

27 Recall that the sum here runs over all m G {l,...,n}, m 7^ I. Likewise, the final product is for all 
k G {1, . . . , n} except I and m. 
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homogeneity zero in all other points. Since the first factor in the sum is everywhere finite, 
the only possible poles in a m come from the final product and so occur when p m collides with 
some other marked point with k,m ^ I. But for any given k, m (say m = 2 and k = 3), 
it is easy to check that the singularity cancels in the sum. Therefore (4.35) has no poles in 
any of the a m (and hence none in any of the cr^.). But by Liouville's theorem, a function 
homogeneous of degree zero that is everywhere holomorphic on a compact Riemann surface 
must be constant. Quite remarkably, we have learnt that (4.35) is completely independent 
of all the marked points except for <7j. Finally, since (4.35) is both a (l,0)-form in 07 of 
homogeneity 2d and is linear in the infinity twistor ( , ), we see that 

\- (Afo), %J) n yr 1 (arl) -j-r (kl) (X(ai),d\(ai)) 

% ; ilm) 11 (a r m) 11 (km) " (a^) ' 1 ' > 

This is exactly E$ in (4.34). 

With this simplification understood, combining (4.31) with (4.33) shows that the picture 
changing operators give 

( n T (Pi) ) = i ^ xd L = det '( HV ) ( 43? ) 

as the factors of (cruder/) cancel. This is exactly the contribution of the conjugate worldsheet 
Hodges matrix in (1.5), again represented by the specific case that we compute the first dxd 
minor (the leading principal minor). 

We now address an issue that may have been puzzling some readers. In the above, we 
implicitly chose to insert the picture changing operators at d of the same points as the matter 
vertex operators. Although this was the minimal choice, was it really necessary? In fact, 
as in usual superstring theory, the picture changing operators may be inserted at completely 
arbitrary points on the worldsheet, and these locations are not integrated over. Rather than 
repeat the standard abstract argument for this (for which see [54, 85]), we shall show directly 
that despite appearances, (4.37) is in fact completely independent of the choice of d insertion 
points. 

As a warm-up, it is easy to see this claim is certainly true when d = 1, for then det(IHI v ) = 
EI]^ = — (A(o"i), dA(<7i))/(crid<7i) and the Vandermonde determinant is unity. Since Z(a) = 
Aa^ + Ba± at MHV level, this becomes simply —{A,B) which is obviously independent of 
the insertion point. 

For the general case, note first that (4.37) is homogeneous of degree zero in each of the 
(T/s. The minor of H v itself can have no poles — it is a polynomial in its entries (4.34), each 
of which are everywhere finite. Thus the only possible singularities in (4.37) come from the 
Vandermonde determinant \a\ ■ ■ ■ a^\ 2 in the denominator. This produces a second order pole 
when any pair of insertion points collide. We shall show that this singularity is cancelled by 
a second order zero from H^ xd . 
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To see this, suppose p\ approaches p 2 with their separation measured by any small pa- 



rameter e that has a first-order zero when they collide. Then for m > 3 we have ^ 
^2m + O(s)- Subtracting rows and columns, the numerator of (4.37) becomes 



1 m 



u dxd 



ey 2 * o( e ) 

[fV (Ml 

^22 (23) 
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0(e) 



as e becomes small. Here, 
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(2,d) 

(2d) 



0(e) 

(2,d) 
(2d) 



T v 
Md 



(4.38) 



However, in section 3.4.2 was defined to 



be the limit of as the two points collide, so 



0(e) automatically. Similarly, 



2H^ 2 + : 



(4.39) 



is by definition e 2 times the second derivative of W\ 2 at <7i = 0*2, plus higher order corrections. 
So H^* = 0(e 2 ). Therefore, as p\ —> P2 we can extract a factor of e from the first row and 
a separate factor of e from the first column in (4.38), showing that the d x d minor of H v 
indeed has a second order zero in this limit. This cancels the second order pole from the 
Vandermonde determinant in the denominator so that (4.37) remains finite. But by the 
permutation symmetry of |H dxd | and \a\ • • -a^\, (4.37) cannot have any poles in any of the 
worldsheet coordinates. Again, a function of degree zero that has is globally holomorphic on 
a compact Riemann surface must be constant, so (T(cri) • • • Y(<7d)) is completely independent 
of the insertion points, as expected for picture changing operators. 

Above we obtained a representation of det'(IHI v ) in which we computed a minor involving 
the same rows and columns. Once again, we can obtain more general representations, in 
which we compute arbitrary minors of H v , by inserting picture changing operators for the 
two flavours (/ii,^) °f antighost at independent locations. That is, we replace 



T(a) 



T(a,a') = 6(^)( Pl ,Z)(a) x 6(fx 2 )(p 2 , Z){a') 



(4.40) 



to compute a minor of HI V from rows and columns corresponding to the independently chosen 
insertion points of 5(^i) and 6 {^2)- As before, since det'(HI v ) is competely permutation 
symmetric in — indeed, completely independent of — all insertion points, there is no real 
difference between any of these cases, although a judicious choice may help simplify some 
calculations. 

The conjugate Hodges matrix appeared to be the most complicated ingredient in the 
gravitational scattering matrix as presented in [32]. Quite remarkably, it has turned out to 
be one of the simplest. 



4.4 The tree-level S-matrix 

Combining the correlation functions (4.25) & (4.37) with the remaining integral over the zero 
modes of the YZ system — i.e., the space of holomorphic maps Z : S — > PT — and dividing 
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by vol(GL(2;C)) to account for the zero modes of the ghosts associated to the worldsheet 
gauge theory, we have found that 





d 4(d+l)|8(d+l) Z l( " f 

— — — — det'(H)det'(H v ) TT / / li (Z(^))(^da 4 ) . (4.41) 




Recalling the genus zero relation d = k + 1 between the degree of the map and the N k MHV 



level, this correlation function is exactly -M n ,k a s defined in (1.5). Summing over all d > 
and allowing all n > 3 yields the complete tree-level S-matrix of N = 8 supergravity. The 
ability to reproduce this formula for the complete classical S-matrix is a highly non-trivial 
test of our claim that the worldsheet model proposed in section 3 does indeed describe N = 8 
supergravity. 

5 Discussion 

We have shown that the worldsheet theory defined by the action (3.32) and BRST opera- 
tor (3.24) provides a twistor string description of M = 8 supergravity. The model depends on 
a choice of infinity twistor, and different choices lead to N = 8 supergravity on flat or curved 
space-times, with the R-symmetry gauged or ungauged. We showed that in the flat space 
limit, g = worldsheet correlation functions in this theory generate the complete classical 
S-matrix of N = 8 supergravity, in the form discovered in [32] and proved to be correct 
in [33]. By interpreting M = 8 supergravity as a twistor string, the present work supplies the 
theoretical framework to explain why this form for gravitational scattering amplitudes exists. 

The ideas presented here suggest many avenues for further exploration. Let us conclude 
by discussing some of these. 

5.1 The SL(2;C) system 

The most immediately important issue is to properly understand the role of the worldsheet 
gauge theory. In the present paper, our primary concern was to reproduce the tree- level 
S-matrix (1.5). At g = 0, the rank 2 bundle C 2 <%> C is uniquely determined by degree of C. 
Consequently, the holomorphic GL(2; C) bundle has no moduli at g = and the antighost M 
has no zero modes. We were able to account for the zero modes of N rather naively, taking 
the quotient of the zero mode path integral by the obvious GL(2; C) action, or 'dividing by 



At higher genus, holomorphic bundles do have a non-trivial moduli space even for a fixed 
curve S, and this moduli space has been extensively investigated from many points of view 
in both the mathematics [76-80] and physics literature [81-83]. In particular, when g > 2 
the moduli space J4 of stable holomorphic SL(2; C) bundles has dimension 3g — 3, while that 
of stable GL(2; C) bundles has dimension (3g — 3) + g with the extra g corresponding to the 



vol(GL(2))'. 
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Picard variety of C A dense open set of N may be identified with the Teichmiiller space 
of S. We have repeatedly mentioned that our twistor string does not involve worldsheet 
gravity and so its path integral does not automatically include an integral over the moduli 
space of Riemann surfaces. Nonetheless, it does know about Teichmiiller space via the moduli 
space of the rank 2 gauge bundle associated to P. In fact, there is even a natural isomorphism 
between TS) and -ff 1 (3\f, TJ$) (see e.g. [76, 79, 83]), so that deformations of the complex 

structure of £ and of the SL(2) bundle are to some extent interchangeable. The mechanism 
by which this is realized in the current context, and the implications for the twistor string, 
cry out for a better understanding. 

A closely related issue is the apparent absence of vertex operators inserted at a fixed 
point p£X, rather than on a fixed section X <—> X. A proper understanding of the SL(2; C) 
system should include an operator which creates a puncture on £ to which our vertex operator 
is attached. Including such operators should amount to allowing (parabolically stable [86]) 
holomorphic SL(2; C) bundles that have simple poles at pi G S, such that the monodromy 
of the associated flat connection is in a fixed conjugacy class Q C SL(2;C). With n such 
punctures, the moduli space of such meromorphic bundles has dimension 3g — 3 + n. 

As an obvious application, vertex operators associated to punctures on X are likely to 
be important if one wishes to have a worldsheet description of factorization [55, 59]. The 
formula (4) for the gravitational scattering amplitudes was shown to obey the expected fac- 
torization properties in [33]. However, the derivation given there was rather involved, because 
by necessity it dealt with the path integral after integrating out everything but the Z zero- 
modes. By working directly with the vertex operators, one should be able to provide a simpler 
proof (following the general pattern in string theory), as the terms that may become singular 
in the factorization limit are isolated more cleanly. 

5.2 Higher genus 

The discussion of section 5.1 has an immediate corollary that perhaps bears some relation to 
the debate about whether M = 8 supergravity could be perturbatively finite [87-93]. Usually, 
string theory is UV finite because the worldsheet theory is modular invariant. We do not 
integrate over Teichmiiller space, but rather over its quotient by the mapping class group. 
This renders harmless any potential divergence as Im(r) — > 0, and this potentially dangerous 
region becomes (real) codimension 2 rather than codimension 1. In the theory studied here 
though, worldsheet gravity is replaced by a worldsheet gauge theory. As mentioned above, 
we expect the path integral to involve an integral over the moduli space of stable SL(2; C) 
bundles, not the moduli space of curves. Although many aspects of the gauge theory (such 
as the symplectic form) are invariant under the mapping class group, it is not clear that 
we should really expect to take the quotient by modular transformations. If not, then the 
potentially dangerous region Im(r) — > is still present. Of course, it is perfectly possible 
that the integrand still happens to have no singularity here — and indeed we should expect 
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this at low genus — but this requires calculation 28 . We are unable to offer the usual string 
theoretic guarantee that there is simply no place for UV divergences to arise. 

Whatever the fate of N = 8 supergravity at higher loops, the current consensus is that we 
do not expect any UV divergences when g < 7 [87-93] . What prospect does the twistor string 
have for computing these 'intermediate' loop amplitudes? Hopefully, the above discussion has 
made it clear that we cannot give a proper answer to this question without first understanding 
the role of the worldsheet gauge theory. Nevertheless it is clear that many properties of these 
amplitudes are correctly reflected by the worldsheet theory. In particular, the zero modes 
of the Pj- and //^-systems will yield higher loop Hodges matrices that have the correct 
dependence on the infinity twistors [ , ] and ( , ) required by factorization (at least in the 
generic case with d > 2 — 2g; in general we would need to work with n 'fixed' vertex operators 
and an appropriate number of picture changing operators for the /37-system). The factors 
of (ij) appearing in these matrices at g = naturally generalize to the appropriate Szego 
kernels at higher genus, while the Vandermonde determinants coming from the correlation 
function of insertions fixing the zero-modes will involve a basis of holomorphic sections of C 
over a genus g curve. All these ingredients can be written in terms of (higher-order) theta 
functions. See [94] for a related discussion in the context of the original twistor string models. 

Even if successful, it is doubtful that the twistor string would reproduce even one-loop 
amplitudes in a form that permits direct comparison with known results in the literature [28] 
(though some of the expressions found in [95-97] may be closer). A direct assault on the 
resulting integrals is unlikely to be successful; the integrals over the moduli space of higher 
degree twistor curves is challenging even at g = [98]. The most promising approach is 
probably to check that the resulting expressions have all the correct factorization properties. 

5.3 Boundary correlation functions in AdS4 

In this paper, we concentrated on taking the flat space limit so as to extract gravitational 
scattering amplitudes and make contact with the known literature. However, the theory is 
equally capable of describing supergravity or gauged supergravity on AdS backgrounds — we 
simply keep the infinity twistor or infinity supertwistor non-degenerate. 

On anti de Sitter space, the natural observables are not scattering amplitudes but rather 
boundary correlation functions. These are obtained by choosing the external wavefunctions 
to represent bulk-to-boundary propagators, that is, solutions of the free equations of motion 
on AdS, with a singularity on the conformal boundary. Such wavefunctions have a very 
simple description in twistor space, known in the twistor literature as 'elementary states' (see 
e.g. [99]). For example, consider the elementary state 



Another intriguing but very speculative idea would be that the theory allows us to take the quotient by 
a g-dependent congruence subgroup of the mapping class group that becomes trivial when g is greater than 
some minimum value g , signalling the onset of UV divergences. 
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representing a scalar field in twistor space. If the line AB is chosen to obey (A, B) = 0, then 
it lies at infinity. In particular, if X is the non-degenerate infinity twistor associated to AdS4, 
then this twistor line represents a point y on the thee dimensional conformal boundary. Using 
the standard incidence relation p a = x aa \ a , the Penrose transform of (5.1) appropriate to 
AdS4 is 

K[r fl) _ , (Z>dZ) _ / (l + A*»)(AdA> 



A-Z B-Z J (A & x«<* + A a )X a (BpxM + Bf>)\ p 

:i + Ax 2 ) 



■ ~-y) 2 

where we used the non-degenerate infinity twistor in the measure (Z, dZ). This is the bulk to 
boundary propagator for a scalar field, written in the coordinates where 

dS " (1 + Ax2)2 V- 6 > 

is the AdS4 metric and where (x — y) 2 is computed using the flat metric. 

Using states such as (5.1), it should be possible to use the formalism of this paper 
to compute arbitrary n-point boundary correlators, again in the form of an integral over 
the moduli space of degree d curves in CP 3 , at least at g = 0. One obvious feature is 
that the (n — d — 2) x (n — d — 2) worldsheet Hodges' matrix and the d x d conjugate 
Hodges' matrix will combine into a single (u — 2) x (n — 2) worldsheet matrix, with the 
off-block-diagonal terms being proportional to the cosmological constant. These terms arise 
because with a non-degenerate infinity twistor, both the pp- and yZ-systems have cross- 
contractions between the Oh vertex operators representing the external states and the picture 
changing operators T. Indeed, one could anticipate this happening. The generalization of 
Hodges' MHV amplitude to the worldsheet Hodges' matrices was deduced [32] starting from 
the observation that factorization requires the n-particle N d_1 MHV flat space tree amplitude 
to contain n — d — 2 powers of [ , ] and d powers of ( , ) when written in twistor space. But 
with a non-degenerate infinity twistor these two objects are really equivalent. 

In the twistor string, as in usual string theory, factorization of scattering amplitudes is 
closely related to collision of vertex operators on the worldsheet [33, 36, 39]. Factorization of 
boundary correlators in AdS has been investigated recently in [100-103], where it is shown 
that (tree-level) Witten diagrams in AdS obey a natural generalization of BCFW recursion. 
An important observation related to this is that the structure 

D 3|8 Z a — A M h {Z Y +tZ n ,...,Z)A Mr(Z, ...,Z n ) (5.4) 

of the BCFW recursion relation in twistor space [37, 38] is completely conformally invariant 
when expressed in terms of external 'twistor eigenstates' 

tH(Z) = t l8 (Z, Zi) = I ^ A ^ 4|8 (Z 4 + sZ) . (5.5) 
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In particular, the infinity twistor arises only via the three-point functions that seed the 
recursion relation. We can obtain AdS boundary correlators by integrating (5.4) against 
appropriate boundary elementary states hi(Zi). Thus, in twistor space, the BCFW recursion 
for Witten diagrams in AdS is exactly the same as BCFW recursion for flat space amplitudes. 
Only the translation back to momentum space (associated to the boundary <9AdS) and the 
three-point functions are different. It would be fascinating to relate these observations to the 
structures of Witten diagrams found in [100-103]. Of course we are limited to the case that 
the bulk AdS space is four-dimensional. 

Finally, via analytic continuation to dS, boundary correlators of gravitational modes on 
AdS may even have cosmological applications [104]. The ideas presented here may provide 
a way to extend the calculations of [104] to higher-point functions. The (n > 3)-point 
gravitational wave power spectrum is admittedly a rather esoteric cosmological observable! 

5.4 Other issues 

We briefly mention various other issues. 

Firstly, the theory we have presented is purely chiral really provides a top holomorphic 
form on the moduli space. It is this 'scattering form' 1 that was found in [32]. To recover the 
actual scattering amplitudes we must still pick a 4d-dimensional 29 real integration cycle on 
which to integrate this form. When d = 1 and g = 0, the moduli space is simply complexified 
space-time and the appropriate integration cycle is just a copy of real Minkowski space. For 
higher degrees the appropriate contour is less easy to define. One possibility, suggested in [12] 
and hardwired into Berkovits' model [56], is to pick real structures 30 T\ on £ and T2 on CP 3 
and ask that the map is equivariant in the sense that Zot\ = T20Z. In ultrahyperbolic space- 
time signature, these real structures fix an S l equator on S at g = and an MP 3 real slice 
of twistor space. However, some care is needed in the interpretation of wavefunctions on real 
twistor space (see e.g. [37] for a discussion). Other integration cycles of interest include those 
that compute factorization channels of amplitudes, ultimately yielding 'leading singularities'. 
It would be good to know whether the twistor string naturally picks a preferred integration 
cycle for us, or whether this is additional data that must be specified. 

In this paper, although we identified the relevant transformations of X that were being 
gauged, we did not attempt to write down a classical action theory that realized this gauge 
symmetry off-shell. Instead, we moved right away to a gauge fixed model together with 
its ghosts and BRST symmetry. It would be interesting to construct the unfixed model, 
particularly is this would likely shed further light on the role of the GL(2)-system. Such a 
model would appear to involve 2 charged gravitinos and 2 charged gauginos in addition to 
the GL(2) gauge fields. 

29 This is in the case that the wavefunctions are represented in terms of Dolbeault cohomology classes 
ff' ' 1 ' (OT, 0(2)). A description in terms of sheaf cohomology would require us to pick a (4d+ n)-dimensional 
cycle. 

30 Recall that a real structure is an antiholormorphic involution squaring to the identity. 
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Next, the vertex operators Oh and Oh that we obtained are natural analogues of Neveu- 
Schwarz sector vertex operators in the superstring. It is important to know what, if anything, 
the Ramond sector could be in the present context. Unlike conformal gravity modes in the 
original twistor string, we would not expect Ramond sector operators to be generated at g = 
if they are not present (pairwise) in the external states. If they exist, their role at g > 1 is 
clearly important to understand. 

We saw in section 4.2.1 that, when restricted to constant maps, the string field theory of 
our model is the twistor action for self-dual N = 8 supergravity found by [66] . The string field 
theory of the full model should thus include a further term representing worldsheet instanton 
contributions. Presumably, only the degree 1 instantons need be included, as is the case 
in usual string theory [105] and as in the analogous twistor action for N = 4 super Yang- 
Mills [3]. In our context, these would represent off-shell gravitational MHV vertices. This 
strongly suggests that despite the difficulties [22] with Risager recursion for gravity [20, 21] 
an MHV formalism for gravity exists. It would clearly be of great interest to find a twistor 
action for non self-dual gravity. The deformed worldsheet action (3.51) perhaps provides a 
good starting-point. See [24] for an earlier attempt to construct a twistor action for gravity. 
An important step in the right direction has recently been taken in [40]. 

Last but not least, it would be very interesting to revisit the potential existence of a 
twistor string for pure N = 4 super Yang-Mills in the light of this paper. One approach might 
be to try to understand the meaning of the duality between colour and kinematics [106] in a 
twistor framework. This duality has certainly lead to great progress in the computation of 
multi-loop gravitational amplitudes in momentum space, typically with u = 4. The similarity 
between the twistor action (4.28) for self-dual gravity and (4.29) for self-dual Yang-Mills is 
surely no coincidence. Yang-Mills amplitudes are completely permutation symmetric in the 
external states provided we include their colour factor. Perhaps they also admit a Hodges 
matrix form. 
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A Some properties of algebraic /37-systems 



In this appendix we will compute some correlation functions of operators in ^7 systems that 
are ingredients in computing the twistor string worldsheet correlator (4.1). Nothing in this 
appendix is new — all (and much more besides) may be found in [107] and in section 10 
of [55], which we follow closely. 

In constructing the twistor string theory, we imposed no reality conditions on the world- 
sheet fields (see the discussion in section 5.4). Thus the path integral over these fields should 
be understood as a formal algebraic operation. This is exactly the usual case for Berezin 
integration of fermionic variables, and so the discussion of [55] proceeds by relating integrals 
over bosonic fields to integrals over fermionic fields that are easier to understand. Thus we 
consider the path integral over anticommuting fields that we call b and c with action 

S bc = ^J^dc. (A.l) 

The result of this path integral depends on the quantum numbers of these fields. Without 
loss of generality, we can assume that 

c€^°(S,L) and b G n°(£, if E ® L" 1 ) (A.2) 

for some line bundle L, and we assume the <9-operator in the action in (A.l) acts appropriately 
on sections of L. Then zero modes of c are globally holomorphic sections of L while zero modes 
of b are globally holomorphic sections of K% (g) L _1 . By Serre duality, this is iT (£,L). In 
the case that L is a spin bundle, so that L 2 = (generically) neither field has zero modes 
and the be path integral is 

j V(b, c) exp (~J bde^j = det(0^ /2 ) , (A.3) 

or in other words the determinant of the Dirac operator on S. As explained in [108, 109] this 
may be written in terms of the Riemann theta functions associated to £ and the choice of 
spin structure. When g = 0, we may take it to be a constant. 

For any other choice of L, at least one of b or c will have zero modes. By the usual rule 
j d6 ■ 1 = of Berezin integration, the path integral (A.3) vanishes. To obtain a non-vanishing 
result, we must insert exactly enough fields to absorb the zero modes. For simplicity, let us 
suppose that c has some number m of zero modes, so that we may expand it as 

m 

c(z) = CjVi(-g) + non zero- modes (A. 4) 

i=i 

where Cj are anticommuting constants and the ^ form a basis of H°(T,,L) (here written in 
terms of a local coordinate z G U C S). With m insertions of c, the path integral becomes 

J V(b, c) c(zi) ■ ■ ■ c(z m ) exp J bdc S j=det'(B L )x J dQ cq{z\) ••• c (z m ) (A. 5) 
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where cq(z) = ^ Ci^i(z), and where the determinant is provided by the path integral over the 
non-zero modes of the ftc-system. Since the c's anticommute, (A. 5) must be antisymmetric 
under the exchange of any pair of insertion points z% and zj. It must also be holomorphic in 
all of these insertion points. Thus we find 

j V(b, c) c{z x ) ■ ■ ■ c(z m ) exp (-^J bdc^j = det'(B L ) x det(Y) , (A.6) 

where Y is the mx in matrix with entries Y, L j = ^i(zj). Equation (A.6) is the standard result 
for fermions. 

In Berezin integration, if c(z) is fermionic then 5(c(z)) = c{z). In addition, because 
e r = 1 + r if t 2 = 0, we can represent S(c(z)) in integral form as 

8(c(z))= I drexp(rc(z)) (A.7) 



where r is an auxiliary anticommuting variable. This is clearly analogous to the usual inte- 
gral representation of the Dirac 5- function. Following [55] we thus introduce m = /i°(S,L) 
such constant anticommuting variables (ti, . . . , r m ) and let b indicate the collection of fields 
(b(z); n, . . . , T m ). We also introduce the extended action 



m 

J s i=l 



and the extended path integral measure 

T>(b, c) = V(b, c) n . z . m dci • • • dc m dri • • • dr m , (A. 9) 

where T>(b, c) n . z . m . is the measure on the infinite dimensional space of non zero-modes. Com- 
bining (A.7)-(A.9) we see that the path integral (A.6) may be rewritten as 

V(b,c) e~ Sbc Jp( c ^)) = / P ^' c ) e_S?c ( A - 10 ) 

8=1 

in terms of the extended set of fields and action. The virtue of thinking about (A.6) in this 
way is that we have changed a path integral with insertions into a simple path integral over 
a Gaussian action. 

It is now straightforward to understand the bosonic case that is actually needed in sec- 
tion 4. Suppose P and 7 are fields on S with exactly the same quantum numbers as b and 
c, except that they are commuting fields. In the case that 1? = ifs, in contrast to (A. 3) we 
have 

m7)exp (__L^5 7 ) = _J_ (A . n) 
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giving the inverse of the determinant, as is familiar from Gaussian integration 31 . When L is a 
more general line bundle such that 7 has zero modes, the path integral diverges (or, without 
a reality condition, is ill-defined) because of the integration over these zero modes. They can 
be fixed by inserting (^-function operators, and again we represent these in integral form as 

8(j(z)) = J dtexp(Hz)) . (A.12) 

Constructing an extended action and path integral measure as before, but now with com- 
muting variables, we again convert multiple insertions of such 5-function operators into a 
Gaussian integral. We thus find 

fvw,,) e*p(-i/^5r) g'M*)) = de , ( gJ det(Y) ^ 

where Y is the same matrix of zero modes as before. Notice that with our formal algebraic 
treatment, there is no modulus sign on the determinants on the right hand side. Notice also 
that if j(z) represents a section of L, then £(7(2)) should transform as a section of L . Both 
sides of (A.13) transform as sections of ®iL~ 1 \ Zi . 

— 1/2 1/2 

In the main text, we will be interested in the cases L = <S> £ and L = ® £, 

where £ is a line bundle of degree d. In particular, when g = 0, C is uniquely determined to 
be CP i(d). The appropriate zero modes are then 



a QLi . . . a °Ld+l 



^) = ~ 1 A U/2 fOT K Z /2 ®£ , s 

(ada) 1 / 2 (A.14) 
y i (a) = a^---a^-i( fT dcj) 1/2 for K+ 1/2 C 

where i runs over all possible choices of the indices a A+l or a 1; . . . , a^-i; respectively. 

Inserting these zero modes into Y in (A.13) gives equation (4.13) for the zero modes of 
each flavour of the worldsheet fields 7 a a, and (4.31) for the zero modes of each copy of the 
worldsheet field fj, a . Recall that det'(5i) is a constant at g = 0. 

Far more information about correlation functions in algebraic (3j systems can be found 
in [55]. 



31 With no reality condition on the /37-system, this is really a definition of what we mean by the Gaussian 
path integral. See sections 3 & 10 of [55]. 
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